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SPECIFICATION 

Electronic Version 1.2.8 
Stylesheet Version 1 .0 

Breast Cancer Risk Analysis and 
Computer- Aided Diagnosis 

Background of Invention 

[0001] FIELD OF INVENTION 

[0002] This invention relates to systems and methods for analyzing digital 

mammograms, and, more particularly, to systems and methods for performing a risk 
analysis and computer-aided diagnosis using digital mammography data. 

[0003] CLAIM OF PRIORITY 

[0004] This application claims priority from U.S. Provisional Patent Application Serial No. 
60/334,643 filed November 30, 2001 . 

[0005] GO VERNMENT SUPPOR T 

[0006] This invention was supported in part by the Department of Health and Human 
Services, National Cancer Institute, National Institutes of Health, under Grant No. 5 
R21 CA79947, and (2) the Naval Medical Research and Development Command funds 
managed by the Henry M. Jackson Foundation for the Advancement of Military 
Medicine contract # 600-0675299-169. 

[0007] BACKGROUND 

[0008] Breast composition, a breast cancer risk factor, is correlated with some of the 
known risk factors and is independent of others. Understanding these interactions 
provides a mechanism for measuring the dynamics of breast cancer risk. The present 
method provides a mechanism for fusing risk analysis with automated computer aided 
detection in digital mammography. 
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[0009] It has been demonstrated that regular screening by mammography reduces the 

mortality rate associated with breast cancer (here in "BC"). However, there are debates 
as to whether the screening benefits are applicable to women of all age groups, and 
as to when regular screening should begin, sometimes argued at 40 years of age, or 
more often 50 years of age, since this is where the mortality reductions have been 
previously observed. However, current research indicates benefits for those that begin 
regular screening at 40 years of age. Setting this argument aside, it is clear that a 
serial electronic database will contain many years of clinically normal images for most 
women, and many years of normal images prior to cancer diagnosis for a small 
minority of women. The magnitude of this data suggests that reliable automated 
methods of serial analysis should be useful. 

[001 0] The present invention relates to the risk factors related to the occurrence of BC, 
examining factors that influence serial tissue change and integrating the associated 
influences observed in the corresponding breast radiographs. The present invention 
advances the art in real time automated analysis in digital mammography (DM) with 
special considerations for serial applications, where risk factors and other patient data 
are fused with the image data for improved decision making. 

[001 1] Full field digital mammography (FFDM) is a newer imaging technology, where the 
digital mammogram is acquired directly in electronic format (film-less) as opposed to 
digitizing screen-film images for processing purposes. This digital system, which is 
intended to replace traditional film-screen mammography, has the potential to 
facilitate various teleradiology applications including electronic image archiving and 
retrieval, remote location analysis, and computer aided diagnosis (CAD) including 
serial analysis. The disclosure provided herein is applicable to both digitized film- 
screen images and direct digital images acquired with FFDM systems. 

[0012] 

The breast radiograph is a rather crude abstraction of the 3D volume. 
Incorporating risk information with image measurements brings more information to 
task when designing CAD methods. The FFMD system allows (as does digitized film- 
screen images) the direct connection of the digital image data with the other personal 
patient data (risk factors) in a form that is easily electronically accessible for real time 
automated applications. Preliminary evidence indicates that in some respects the 



APP ID=1 0065929 



Page 2 of 91 



images are similar with digitized images acquired from screen-film in that cancer 
detection rates are similar. 

[001 3] A short review of Wolfe's seminal work is helpful. The assertion was put forth that 
it may be possible to correlate breast patterns with the predilection of BC (this is 
Wolfe's driving hypothesis). Early on Wolfe observed a relationship between the 
prominence of ductal patterns and BC in breast radiographs. The associations with 
lobular tissue (Wolfe referred to this as alveolar tissue), ductal prominence and risk 
were investigated. Wolfe believed that fat, connective tissue and the epithelial 
elements are seen as dysplasia (increased radiographic density constituted by 
interlobular connective tissue) on a mammogram and that periductal fibrosis is 
represented by prominent ductal patterns. Dysplasia is now an obsolete term, which 
was used in the past in reference to fibrocystic disease of the breast; for historical 
reference we use this terminology if the respective authors also used it, which now 
can be considered in general as increased density. Wolfe made the observation that 
for normal patients the ductal appearance becomes more prominent in size and 
frequency with age and decreases with increasing parity, while lobular tissue 
decreases with age. Similar trends for patients with benign disease were also noted. 
However, in women with malignant disease, the ductal patterns are more prominent at 
an earlier age, in addition to having higher proportions of lobular tissue (dense 
tissue). The correlation of changing tissue patterns with age that Wolfe observed did 
not take into account the fact that his early images were film mammograms that had 
higher contrast than the subsequent images that were made with xeroradiography 
which had a wider exposure latitude that can make dense tissues fade. 

[0014] Wolfe made the observation that in some cases, even in the absence of a 
recognizable mass, cancer diagnosis (especially the scirrhous variety with 
predominantly connective tissue character) might occasionally be made by the 
appearance of a unilateral group of dilated ducts. In the past the presence of 
asymmetric breast tissue, especially if increasing in size has been regarded as a 
mammographic sign of malignancy. The American College of Radiology Breast 
Imaging Reporting and Data System (BI-RADS) lexicon has defined asymmetric breast 
tissue relative to the contralateral breast as " a greater volume of breast tissue, 
greater density of breast tissue, or more prominent ducts. Asymmetric breast tissue 
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has been reported to be depicted on 3% of mammograms and is nearly always benign. 
Therefore, a stable finding that fulfills the BI-RADS definition seldom warrants a 
biopsy. Sometimes focal asymmetric densities may represent masses with borders 
that are ill defined or are obscured by surrounding fibroglandular tissue rather than 
asymmetric tissue. These need additional evaluation especially if palpable. 

[001 5] In subsequent work the idea of increased cancer risk and ductal pattern was 

expanded. Wolfe notes that the prominent ductal presence is an abnormality caused 
by sub-epithelial collagen deposits and that the cross-sectional diameter of the ducts 
is directly proportional to the amount of collagenosis. The increased involvement and 
severity have a very similar age dependency as does breast cancer. A detailed 
description of ductal image appearance was also provided. Briefly, these are linear 
bands of increased density that fan out from the subareolar region, and have a 
discontinuous winding or crossing pattern in a severely involved breast. The presence 
of cancer alters the duct pattern distribution, resulting in enlarged and winding ducts. 

[001 6] Following this early work, Wolfe developed a detailed classification scheme 

consisting of four breast patterns that are correlated (this is the driving premise) with 
BC risk. In this scheme the breast composition is classified by appearance (normally 
classified by a radiologist's assessment). In ascending relative risk order, the 
designations are Nl , PI , P2, and DY with disregard to abnormalities, vascularity, or 
any other associated characteristics or risk factors. Here the relative risk means as 
compared with Nl . The Nl , normal or negative, is an image comprised primarily of fat 
with no visible ducts (lowest risk). PI , for prominent duct pattern, low risk, consists 
primarily of fat with visible ducts occupying up to one-fourth of the volume in the 
subareolar region (low risk). These ducts appear as cord like or beaded linear 
structures extending into one quadrant in some cases. P2, the prominent duct 
pattern, occupies more than one fourth of the volume of the breast and may involve 
all of the parenchyma in its severe form. DY, the dysplastic pattern, is the fourth 
category, which poses the highest risk, shows increased density greater than that of 
fat with absence of a prominent ductal pattern. In this study evidence was presented 
that the cancer incidence risk progresses with pattern in a significant manner when 
examining a referral population. It was also noted that not all carcinomas are 
accompanied by prominent duct patterns suggesting that there may be two types of 
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different age related pathogeneses that indicate risk. 

[001 7] In addition to this grouping, a fifth category, QDY, was added for younger women 
below the age of 40-45 to account for physiological changes in the breast due to 
aging. The QDY classification is for women whose radiographically increased density 
is not as severe; these patterns are the portion of the DY images that are most likely 
to regress to a low risk pattern with age. In the same work, two studies were 
presented that corroborate the earlier work in showing a significantly increased cancer 
incidence for women with high-risk patterns. 

[001 8] Another important observation stemming from this work is related to the 

changing role of the high-risk patterns with age. Women with the DY patterns account 
for a large portion of the cancers in the under 50 age group whereas after this age the 
proportion of breasts with DY pattern is significantly reduced, and women with the P2 
patterns account for the larger proportion of cancers, again suggesting the possibility 
of different agents at play. The age related pattern regression is discussed in more 
detail in Part 2. The evidence is almost unequivocal in that the high-risk patterns 
generally regress in time. 

[001 9] In subsequent work, the influence of the aging process on patterns was 

investigated further by studying a referral population. It was found that women 
categorized with either low risk or very high risk patterns at initial assessment do not 
demonstrate much change over time, and that women designated with the QDY 
pattern sometimes regress to a low risk pattern with age. It was also shown that the 
patterns change very little for women after the age of 50. In summary it was found the 
N 1 changes to PI in 1 4% of the cases, PI to N 1 in 7% of the cases, DY regressed to P2 
in 33% of the cases and to PI in 1 1% of the cases, and to Nl in 8% of the cases, and 
the P2 pattern rarely reverted to another designation. The processes of aging and 
involution are discussed again with more detail in Part 2. 

[0020] | n | ater wor | <) Wolfe discusses the elevated incidence risk observed in the P2 and 
DY patterns and that it is extremely unusual to have a cancer develop in a breast once 
classified as Nl . Wolfe also discusses the time interval, on average about 48 months, 
between the negative mammogram and histological proof of carcinoma which 
reinforces the idea that the cancers were not missed at initial screening (this is often 
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termed the masking effect and is discussed in detail below). This work also suggests 
that true DY (not QDYs) patterns rarely change, and if so revert to P2 patterns; thus 
the risk does not decrease with time, since the P2 patterns are stable in time. 

[0021] This early work by Wolfe suggests a significant relationship between the breast 

composition and risk. Additional support was subsequently provided by Wellings and 
Wolfe by correlating histological data with breast composition data from 143 biopsies; 
the highest grades of precancerous epithelial abnormality occurred most often in the 
combined P2DY group and infrequently in the PI category and very rarely with the Nl 
composition. This should be taken with caution, since biopsy material was 
investigated, which is a priori abnormal tissue. The results may not, in general, be 
extended to include normal tissue related patterns without further investigation. 

[0022] There are other methods used to characterize mammographic densities in 

conjunction with or independently of Wolfe patterns. One more currently used method 
involves measuring the amount of radiographic densities in the image; that is 
measuring the tissue that is not associated with fat (radiolucent areas). Another 
method is the BI-RADS system, but there is little work in the area of relating BI-RADS 
image classification with risk. For clarity and where appropriate, these methodologies 
and associated outcomes are addressed separately. For clarity, the Wolfe-pattern risk 
and density-risk are discussed separately when necessary. 

[0023] 

The work presented by Wolfe suggested (a) a measurable risk factor derived from 
the image, (b) the interval between screening sessions may be adjusted according to 
breast composition (risk), and (c) screening practice, in general, may be based on 
parenchymal pattern. This created much debate and early on many attempts were 
made to duplicate and verify Wolfe's pattern-risk premise. The results of this 
verification work provide a less than complete picture, since some studies 
corroborated Wolfe's premise to differing degrees, and yet other research could not 
demonstrate the association between dense breasts and increased cancer risk to the 
same degree or not at all. However subsequent analysis, that reexamined some of the 
earlier pattern-risk work by imposing strict methodological standards, concluded that 
there is an increased risk associated with dense breasts for women of all ages. 
Likewise, additional retrospective considerations indicate that carefully conducted 
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epidemiological studies support Wolfe's original assertions. Even though much of this 
work supports the premise that patterns are related to risk, not all agree that 
screening can be based on pattern (risk) assessment. 

[0024] Much of the research supporting Wolfe's original work does so in varying degrees. 
It is instructive to give a summary of these variations by experimental design. In 
general there are three classes of studies (a) case-control (CC), (b) cohort (C) and (c) 
prevalence (P). In the CC study, cases are matched to controls by some criteria (maybe 
randomly) and compared, which provides Odds Ratios (ORs). Under certain conditions 
ORs approximate relative risk (RR) estimates derived from cohort studies of the same 
population. Incidence (cohort) studies allow for true RR estimates, and prevalence 
studies provide a cross-sectional comparison of those with the disease compared to 
those without the disease at a particular snapshot in time. CC studies also provide a 
relative cross-sectional view. First, the pattern-risk assessment is provided. 

[0025] The earlier prevalence studies (or cross-sectional studies) corroborating the 

pattern-risk assertion show varying degrees of support when stratified by age and or 
pattern. A brief synopsis of this work is given in order to provide a full understanding 
of the variation. 

[0026] Tabar et al make the important observation that the cancer prevalence rate curve 
for the Nl PI (that is the combined Nl + P1 ) group for older women resembles that of 
the younger P2DY (P2 + DY) women with an approximate (displacement to the right) 1 5 
year delay. Other reports also show a diminishing risk with increasing age for women 
with DY patterns. Elsewhere, it has been shown that the P2 pattern may be the highest 
risk category. Other reports provide evidence, to the contrary, that does not support 
the pattern-risk premise for prevalence cancers. 

[0027] Some reports show increased risk for women of all ages with DY patterns and 

increasing risk with increasing pattern. And yet other research indicates a diminishing 
risk with age for dense patterns, with one study showing an opposite trend. Some 
research indicates that the P2 pattern may be the highest risk category. In contrast 
other evidence indicates an increasing risk for the DY patterns with only a slightly or 
no elevated risk for P2 patterns, respectively. In comparison, current research 
involving Singaporean women indicates a high risk associated with a pattern similar to 
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P2 but little risk associated with a classification with DY similar patterns. Other 
research indicates in general that women in the P2DY group are at an elevated risk 
compared with the low risk Nl PI group with significantly increased risk associated 
with the DY pattern compared with the Nl , which is in agreement with subsequent 
work. Evidence suggests that the relationship between high-risk patterns and risk is 
strengthened when accounting for weight and height. Subsequent to the earlier work, 
research by Wolfe and other researchers provide general support for the risk-pattern 
association. One report deserves special attention because it shows an elevated risk 
for P2DY patterns highest for women in the 45-60 years of age; a study important 
because the patterns were recorded four years prior to diagnosis. Subsequent reports 
show slightly elevated BC rates for P2 and DY patterns detected at screening and 
significantly elevated risks for interval detected cancers within 1 8 months of a 
negative screen, and yet other research indicates the pattern-risk assertion holds for 
postmenopausal women only. Bucchi et al find that women with benign breast disease 
are more likely to have high risk Wolfe patterns compared with normal women, and 
have pattern distributions similar to that of women with proliferative and malignant 
disease. 

[0028] Some reports do not support the pattern risk premise. Other work shows 

increased rates for women with dense patterns between 40-54 years of age with the 
risk falling to unity above this age. 

[0029] Gravelle et al provide strong evidence for the pattern-risk association with about a 
4-fold risk increase for women with P2 and DY patterns. One report demonstrates that 
when assessing the incidence odds with patterns taken initially, the BC risks are only 
marginally elevated for the P2 and DY patterns at the end of 8-year screening project. 
However, when accounting for the longitudinally shifting patterns (assessing the 
pattern in the screening round preceding diagnosis) the incidence risks are 
considerably elevated for the high risk group. Although it should be taken into 
account that this group work tracked women initially between the ages of 40-47 years 
of age which is the age group that is most likely to experience serial pattern change. 

[0030] 

The incidence pattern-risk research also shows considerable variability in risk and 
associated follow up time. Krook et al show significantly elevated risk for women over 
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45 for the P2DY group. Research also shows significantly elevated incidence rates for 
the DY class in a four-year (max) follow-up with only slightly elevated P2 odds 
associated with older women. On the contrary, other evidence shows elevated 
incidence rates for both P2 and DY patterns, particularly with the P2 group over a five 
year interval, de Stavola et al (a hybrid study) show an elevated risk (2 fold) for both 
pre and postmenopausal women with P2 or DY patterns at the 5 year censoring time, 
which is in agreement with other research showing increased risk (monotonic increase 
with pattern) with pattern at the five year follow-up. It is suggested that masking may 
be in effect since the proportion of stage 1 tumors detected in the P2DY women was 
less than that of the low risk Nl PI group; this work also shows that by not adjusting 
for age the incidence risk calculation will be underestimated. 

[0031] Briefly, the masking effect implies that cancers are sometimes missed in the more 
dense patterns at initial screening. These missed cancers will then manifest in the 
near future, which will inflate the incidence risk estimates in the short term. Work by 
Egan shows an elevated incidence for dense patterns at 36 months with the 
conclusion that the risk is an artifact of masking. Egan and McSweeney show an 
elevated risk in the P2DY patterns for up to 48 months, which then falls to unity. This 
work is interesting in that it gives lead times (time between normal and abnormal 
reading) for the various patterns, as 43,44,33, and 29 months, for the Nl , PI , P2, and 
DY patterns, respectively, and indicates that most (68%) of the developing neoplasms 
in P2DY breasts occur by the age of 55. Other work that provides little support for the 
pattern-risk premise include, or for five through nine-year follow-up periods. 

[0032] Moreover, it is known that the BC incidence rates vary geographically. Asian 

women are at a lower risk compared with women from the west. Evidence indicates 
that women from India have a higher incidence of low risk Wolfe patterns compared 
with Caucasian women, which provides an indirect link with pattern and risk. However, 
this should be taken with caution. Work based on applying BIRADS tissue 
compositions indicates Asian women are more likely to have BIRADS 3-4 compositions 
compared with Caucasians while Native Americans are more likely to have less dense 
patterns and much lower BC rates compared with Caucasian women. It is the former 
relation that does not fit well with the incidence rates when considering density as a 
risk. Although, other research indicates that Japanese women are less than half as 
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likely to have P2 patterns and 4 times more likely to have N1 patterns compared with 
British women but both groups have similar DY proportions. However, these 
discrepancies may be due to both the tissue rating system and body morphology. In 
support of this premise, other research based on quantitative density measurements 
shows that the mean dense area for Asian women is less than that of Caucasian 
women, but because of the relatively smaller breast size of Asian women, the 
associated percent density is larger thus indicating the confounding influence of 
breast size and density. Research in breast carcinoma among immigrants in Australia 
and Canada shows a shift towards the rate of the destination country, resulting in the 
conclusion that environmental and lifestyle factors associated with the new place of 
residence influenced the incidence of breast carcinoma. An interesting observation 
was that the rate of incidence and mortality among immigrants took place not only 
between from low to high-risk countries but also vice versa. 

[0033] Some researchers have found that relating the amount of radiographic density in 
the image to risk (density-risk assessment) is a more productive methodology. The 
loose definition for radiographic density is any image area that is not fat or not 
radiolucent. Roebuck et al provide a more complete definition for radiographic 
density. It may be the case that density measurements are a more reliable indicator of 
risk. The work discussed in this section, for the most part, is cross-section (case 
control) in nature unless otherwise indicated. 

[0034] Early work by Boyd et al showed that women with at least 1 0% dense tissue which 
he termed dysplastic tissue, were at an elevated risk, and women with at least 75% 
dense tissue were at significantly increased risk; although, this work shows a weak 
connection with ductal prominence and risk. The work by Brisson et al is presented in 
terms of more common radiological terminology when describing densities; in 
conjunction with applying Wolfe patterns, densities are further classified as nodular 
and homogeneous (confluent). Brisson et al show that women with P2 patterns have 
extensive nodular densities while women with DY patterns have extensive 
homogeneous densities as well as nodular densities, where the nodular densities are 
smaller and more concentrated in DY patterns. The relation or correlation with density 
proportions and Wolfe pattern classification is discussed in more detail below. This 
work indicates some interesting points (a) for women under 60 years of age an 
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increase in the percentage of the breast showing nodular densities is correlated with 
an increase in BC risk with a similar, but weaker, trend for homogeneous densities, 
and (b) and increase in the average size or degree of concentration of nodular 
densities is associated with an increased BC risk. The relation of BC risk with nodular 
density percentage is strengthened when weight and height are taken in account, and 
the risk associated with the percentage of the breast showing homogenous densities 
is not observable until accounting for weight and height. Supporting work shows that 
the BC risk rises steadily with increasing density percentage. Safltlas et al show that 
the risk gradient is approximately linear with respect to increasing densities. Byrne et 
al also show an elevated BC odds with increasing density when adjusting for weight. 
The evidence indicates the risk holds when considering age. Although, one report 
shows elevated risks for women with 25% or greater radiographic densities with no 
trend of increased risk with density after this threshold. Byrne et al also found that 
with increasing absolute area of the breast with dense tissue increased also shows an 
increased risk; compared to women with no breast density, women with a total density 
in cm2: 1 -1 3.9, 1 4-22.9, 23-33.9, 34-52.9, and 53 or more had odds ratios of 1 .48, 
1 .99, 2.08, 3.24 and 3.35, respectively, which is in agreement with other work. The 
elevated risks with increasing total density are similar (central values) to those found 
with density proportions. 

[003 5] Few studies have reported on BC risk and percent density as related to time. One 
report also suggests increasing density gives significantly elevated risks for 1 0 years 
of follow-up. One report using a two class percentage density categorization with a 
cut-off at 2 5% density shows an initial elevated risk for the more dense group that 
increased even further at the 3-4 year interval and then decreased with time, which 
supports the masking premise (or possibly accelerated tumor growth); although, the 
masking effect can not explain the initial elevated risk. Brisson et al using a three 
class density categorization of atrophic (fat), intermediate, and glandular or 
homogeneously dense, categories show a stable elevated risk for women of glandular 
and intermediate breast categories of ages 3 5-49 years at the time of entry. In women 
who were 50-74 years of age at time of entry, an elevated incidence of breast cancer 
was observed that diminished with time. A commentary on mammographic densities 
and related factors is also provided by van Gils. 
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[0036] Some researchers have considered risk, Wolfe classification, and density 

simultaneously. Byrne et al show that within the P2 and DY patterns, increasing 
percent densities corresponds with an increased cancer risk referenced to the Nl 
pattern this research indicates that the raised risk with the high risk the Risk within 
high risk. Moreover, the work by Byrne et al indicates that within a particular high risk 
Wofle pattern increasing density accounts for the increased risk. Other work indicates 
the BC risk is elevated with increasing densities for women with P2 patterns, but 
shows no elevation with DY patterns with less than 45% densities; however, an 
increase in risk with density rising above this threshold is observed. The percent 
density risk relation is further strengthened by research that correlates abnormal 
histological change with increasing percent density that may indicate an association 
with high-risk histological changes. 

[0037] Wolfe Pattern and Density Comparison: only two studies were found that quantify 
the association of Wolfe pattern with density proportion. Brisson et al split the 
association (mean values) into percentages of nodular and homogenous densities per 
Wolfe pattern class: Nl , 0% nodular and .4% homogenous; PI , 13% nodular and 1 .9 % 
homogenous; P2, 45.9 % nodular and 3.4% homogeneous; and DY, 24% nodular and 
36.7% homogenous. Wolfe et al found that Nl breasts contained less than 0.01 % 
densities, 80% of the PI breasts contained less than 25% densities indicating some 
difference in correlation. The work indicates a good degree of correlation for the 
high-risk patterns in that 98% percent of the P2 breasts contained at least 25% 
densities, and all DY breasts contained more than 25 % densities. 

[0038] Some research shows BC risk is not related to breast size with or without tissue 
considerations, or related to total breast area, which is in agreement with other work 
when considering weight and height. Although to the contrary, work by Kato et al 
shows a positive association with BC risk and breast area (or volume) for pre- 
menopausal women, but not for postmenopausal women, where a negative relation 
(minor) holds while other work shows an inverse relation with risk and size that 
disappears when considering composition. Yet, research shows an inverse relation 
with increasing breast size with either high-risk patterns or with increasing density 
with one report showing an exception for pre-menopausal women with P2 patterns, 
and that the relation holds for all age groups. Further analysis, when controlling for 
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density, suggests large breasts are normally more fatty and should pose no excess 
risk, but when large breasts are not fat there is a considerable BC risk. Salminen et al 
also found that the incidence of favorable serial pattern change is influenced by 
weight, but the effect is not significant when considering combined influence with 
weight (BMI), age, and number of pregnancies. In general, evidence suggests that 
there is little association with BC risk and breast size when pattern is not considered 
or association between adult height and BC risk. 

[0039] Related research shows that the total area of the breast strongly correlates with 
the total non-dense area and much less to the total dense area in pre-menopausal 
women. The evidence shows that height is associated with total breast area 
(negatively), area of non-dense tissue (negatively), and percent density (positively) 
after controlling for weight and other risk factors. These researchers suggest that this 
correlation work indicates fat and dense areas should be considered separately in 
etiological studies. Although, this dense non-dense correlation analysis results from 
studying pre-menopausal women, it may be safe to assume it applies to 
postmenopausal women also, since menopausal status is not related to breast size. 
Moreover, other work indicates that breast size is independent of age, parity, age at 
first birth, or hormone replacement therapy. 

[0040] One factor in the incidence based risk estimates is that density proportions (or 

Wolfe patterns) change with time, but it is not clear if the BC risk changes accordingly. 
One report shows that for postmenopausal women that initially have 5-25% densities 
that regress to less than 5%, the risk also diminishes. However, no clear conclusions 
can be made as to women with >25% density at the beginning that regressed. 
Similarly, with the cross-sectional analyses the composition assessment time is an 
important consideration. 

[0041] 

The breast is a heterogeneous composition of adipose tissue, epithelial cells 
(parenchymal), and fibrous connective tissue (stromal), and most breast cancers arise 
from the ductal epithelial cells; since fat is translucent, the density is a result of 
fibrous and epithelial cells. The total density, as opposed to the proportion, that is the 
true measure of risk. The differing correlations of dense and non-dense tissue with 
total breast size aids in quantifying the inverse relation of breast size with high-risk 
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pattern, and also provides support for the total-density hypothesis. 

[0042] The hypothesis that certain tissue types (total) are the foci for cancer development 
is supported, circumstantially, by work that shows women who have undergone 
surgical breast reduction, due to breast hypertrophy, have decreasing risk with 
increasing tissue removal. This association should be taken carefully since the risk 
relation may be related to the condition. Moreover, subsequent research indicates a 
significantly reduced BC risk with breast reduction surgery and that the risk reduction 
is inversely related to the amount of tissue removed, in particular for women over 40 
years of age. The total-density hypothesis is also supported by the work that 
investigated total density and risk. Moreover, the research showing that Japanese 
women (low risk) have less total density compared with Caucasians (higher risk), while 
their associated dense proportions are greater renders additional support. 

[0043] It has been shown that magnetic resonance imaging (MRI) tissue parameters such 
as water content and transverse relaxation time, denoted as T1 , are useful for 
separating low risk from high risk Wolfe patterns, high risk patterns show more water 
content and longer Tl relaxation. Likewise, other work in quantifying the association 
with mammographic density and MR breast imaging parameters also indicates a 
positive correlation with percent density and water content, and an inverse association 
with percent density and T2 relaxation; that is the T2 is longer for fat tissue compared 
with that of fibroglandular tissue. This inverse relation is similar to that found with T2 
relaxati on and the DY Wolfe pattern. Moreover, other work indicates a positive 
correlation between density volume proportions obtained from MRI compared with 
planar projection density proportions measured from standard mammography but 
with marked differences in tissue distributions particularly for dense breasts. This 
work indicates that mammography is limited in its ability to accurately assess breast 
tissue proportions due to differences with the three dimension volume and 
associated, reduced in dimension, x-ray projection. 

[0044] 

The totality of the work discussed above indicates that density classification is an 
important factor to take into account when analyzing mammographic images; the 
proportion of radiographic densities in the image is related to the elevated risk. Some 
researchers have acknowledged the importance of automated density measurements. 
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Moreover, research shows that the quantitative evaluation of densities by planimetry 
correlates well with the radiologist's assessment, Similarly, Boyd et al show that, by 
either the radiologist's estimation or with a computer assisted manual thresholding, 
more proportions of radiographically dense tissue are associated with a higher risk. 
This implies that well designed automated methods have the potential for similar 
performance. 

[0045] 

Magnin et al using texture measures and local statistics classify images into four 
distinct categories of risk with some success and conclude that acontinuous labeling 
may be a better approach. Caldwell et al apply a fractal measure for image pattern 
classification and risk assessment resulting in good agreement with the radiologist. 
Kallergi et al apply a local thresholding technique and region seed growing method 
(two different methods) with preliminary results for Nl and DY classification. Taylor et 
el investigated different measures for estimating texture. This work shows that 
skewness followed by fractal characterization allows the best discrimination between 
dense and fatty type images among the various metrics investigated, although the 
skewness distributions of fatty and dense breasts show some degree of overlap. 
Tahoces et al investigated various metrics including Fourier measures, local-contrast, 
and gray level distribution analysis for Wolfe pattern classification based on the 
analysis of a representative region of interest that resulted in reasonable image 
classification. Byng et al found that combining a global fractal measure with a global 
skewness figure of merit, based on averaging local measurements, results in better 
parenchymal density detection as opposed to either figure alone; both measures were 
found to be correlated with the radiologist assessment. In subsequent work, this 
method was applied for assessing risk that showed both parameters were significantly 
related to risk; although, the correlations with the automated density analysis 
compared with the reader's assessment was less than in previous study. In addition, 
earlier work by Byng et al applied a computer assisted thresholding technique for 
pattern analysis. Karssemeijer applies a KNN classifier based on features found by 
transforming the data into regions with equal distances to the breast skin line. Huo et 
al show that women with dense breasts who have a genetic BC predisposition that are 
at a high risk in accord with the Gail model can be detected with texture measures. 
Heine and Velthuizen show that a rigorous parametric statistical analysis allows the 
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continuous labeling of dense and non-dense tissue following a pre-whitening 
technique applied to the raw image based on 1 / f noise modeling. Sivaramakrishna et 
al transform the data into a normalized local variance representation for automated 
density detection based on thresholding. 

[0046] There is considerable research in investigating the associations of breast tissue 
patterns (and or density) with the other known BC risk factors. This work may be 
important for understanding biological mechanisms or hormonal interactions related 
to BC. The thought behind much of this work is that other risk factors might possibly 
mediate through breast tissue or that dense tissue is a surrogate for other risk 
factors. The idea to keep in mind is that if a particular event or action increases the BC 
risk, then a corresponding increase in density or upward shift in Wolfe pattern may 
follow. On the other hand, if some interaction is assumed to reduce risk there may be 
a corresponding density decrease or downward shift in Wolfe pattern. Indeed, this 
dynamic appears to be at least partially correct, as discussed below and in Part 2, 
which may provide a measurable metric for assessing dynamic risk behavior. 

[0047] The present invention accounts for risk factors that act independently (or not 
correlated to a high degree) of the tissue characteristics as well as those that are 
partially correlated with tissue. These factors will be used to incorporate risk analysis 
into the automated detection methods. From an information theory vantage point, 
independent factors bring more to the decision making power. Since the density 
measurement follows from the image analysis, the identification of other factors that 
are not correlated with density is of prime importance to the analysis. 

[0048] Often univariate and or multivariate methods are employed to study risk factor 
iterations. The former assess the relationship of the event (tissue characteristic) with 
the risk (the variate) singularly. The latter is a stronger approach in the sense that the 
interactions between the variates are also assessed in relation to the event. 

[0049] Evidence indicates that nulliparity (positively), parity, or multiparity (negatively) 
and late age of first child birth (positively) are associated with the incidence of high 
risk breast patterns. Many of the reports indicate that with increasing parity there is 
an associated decrease in the likelihood of high-risk patterns. For clarity, it is 
convention to refer to the Nl and PI patterns as low risk and the P2 or DY patterns as 
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high risk regardless of the experimental outcome. This convention is followed here 
where appropriate. 

[0050] Many of the univariate studies (or similar type analyses) indicate, to some degree, 
that later age of first birth is related to higher prevalence of high-risk patterns with 
some reports showing a weak trend. Early age at first birth is an important 
determinant for the development of favorable breast patterns and this influence may 
remain life-long. Bergkvist et al noted an increase from 1 0.8 % in the proportion of 
women with P2 or DY patterns with age of first birth at 21 compared with 1 8.7 % at 30 
years of age at first birth and 31 .3 % among nulliparous women. These effects were 
marked among oldest women. However, other multivariate work indicates that the 
association with age of first birth is not significant when controlling for parity. 

[0051 ] Similarly, multivariate analyses show that later age at first birth is associated with 
high-risk patterns, and parity is related to less dense or low risk patterns with the two 
factors having independent influences with parity having the greater effect. Bartow et 
al observed a significant association between radiolucent patterns and parous states 
with OR of 2.75 for parous women. Similarly, Gram et al found nulliparous women are 
more than 2 times more likley to have high-risk patterns compared with women who 
gave first birth in their teens, and determined that a multipara women with > 4 
children are 90 % less likely than a nullipara women to have high-risk patterns. It is 
interesting to note evidence suggests the magnitude in shift from high to low risk 
pattern from nulliparity to parity is of the same magnitude as the shift from one 
pregnancy to the next. Gravelle et al found, among nulliparous women , that only 1 /3 
displayed the Nl pattern whereas the Nl pattern was three times more likely among 
parous women. Similarly, de Stavola et al show that parous women are about 30-35% 
less likely to have P2 or DY patterns respectively, relative to nulliparous women, which 
holds for both pre and post menopausal women. In addition to showing that parity is 
related to low risk patterns, regression analysis suggests that the shift from high to 
low risk pattern due to parity is primarily due to the loss of P2 patterns with the shift 
to PI or Nl without preference. This work also suggests that the P2 patterns decrease 
with increasing pregnancies with no further effect on the DY patterns when 
considering 1-5 pregnancies. 
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[0052] 



The connections with density proportions are somewhat similar. Brisson et al show 



some interesting associations (a) increasing number of full term pregnancies is 
associated with a decrease in nodular densities with little variation due to menopausal 
status, (b) a similar trend in the mean percentage of breast showing homogeneous 
densities was shown for pre-menopausal women only, and (c) no trend was apparent 
for the age of first birth for either density and some association with family history 
and density (both types) was indicated for pre-menopausal women. Using multivariate 
analysis Saftlas et al show that percentage of densities decreases with increasing 
parity. Boyd et al compared normal women with 25% densities with women that have 
75% or more densities and found that women with increased radiographic density are 
more often nulliparous. One report indicates that the combined effects of nulliparity 
and density act in synergism indicating that density is more than a surrogate for 
nulliparous status. This work shows parous women with 5-25 %, and >25% density are 
at about a 2.7 and 3.6 RR of BC, respectively, and nulliparous women with >5% 
density have about a 7 fold increase of BC relative to parous women with <5% density. 

[0053] Menstrual Factors: An association with the age at menopause was also made by 

Gram et al indicating that menopause between the ages of 45-52 was associated with 
a 30 % likelihood of high-risk patterns, and a 50 % likelihood of high-risk patterns if 
age of menopause was 52 years or older compared with menopause below 45 years of 
age. 

[0054] Generally, the evidence relating age of menarche with pattern is ambiguous. Gram 
et al show that increasing age of menarche (1 6 years or more) results in double the 
likelihood of having high-risk patterns for pre-menopausal women, whereas an 
inverse relationship was found for postmenopausal women who were 80% more likely 
to have a less dense pattern if their age at menarche was later (past 1 6 years). 
Similarly, later age at menarche was found to be associated with dense Wolfe patterns 
in univariate analysis, which to the contrary only held slightly for postmenopausal 
women in a multivariate work with adiposity explaining the effect in pre-menopausal 
women. The relationship of early age of menarche and low risk pattern is supported 
by earlier work, although menopausal status was not considered. Yet, other work finds 
no connection with age of menarche with pattern. Other work provides a clearer 
interpretation of menstrual history by using the variable years of menstruation and 
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found that the odds of breast cancer increased steadily with the number of years of 
menstruation. Estimates indicate the odds for developing breast cancer nearly 
doubled, somewhat linearly, when the number of years of menstruation increased 
from <2 5 to >40 years. 

[0055] Reproductive factors such as parity, age of menarche, age of menopause, and 
length of the reproductive years have been related to BC risk in that, later age of 
menarche, early age of menopause, either natural or surgically induced, and parity 
provide protective effects against BC. A significant inverse relationship with age at 
menarche and BC risk in pre-menopausal women has been noted with decreased BC 
risk associated with a menarche age of 1 5 years or more compared to menarche at 
age 1 2 and 1 3, respectively; the reported reduced risk indicates ORs of about .45 for 
nulliparous pre-menopausal women and .72 for pre-menopausal women (without 
regard to parity) with a weak risk reduction for postmenopausal women (without 
regard to parity). This does not appear to fit well when compared with the age of 
menarche tissue relations discussed above. To the contrary, evidence (without regard 
to pattern) shows increased BC risk with later age at menopause with an OR of about 
1 .91 among nulliparous women with menopause age of 53 and greater compared with 
those of menopausal age of < 45 years, which does agree with the age of menopause 
tissue relations discussed above. 

[0056] The parity-risk relation may help explain some of the conflicting reproductive- 
pattern associations. Research, comparing uniparous and nulliparous women, 
suggests that early childbirth provides a protective effect for all periods after birth, 
whereas late age of first birth induces a higher risk immediately after birth that 
declines with time. 

[0057] 

The associations of height, weight, and body build with pattern risk analysis 
operate in opposition since tall thin women are more apt to have dense patterns or 
more percentage density, whereas short heavier women are more apt to have low risk 
patterns. Two reports demonstrate that the BC risk associated with weight is not 
apparent unless parenchymal pattern is taken in account, and when considered height 
is no longer important. Brisson et al found that the proportion of women with P2DY 
patterns was 93.6 % among tall and thin women (height > 1 65cm, weight < 55 kg) 
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and 1 8.8 % among shorter and heavier women (height < 1 55cm, weight > 75kg)(57). 
This is in agreement with subsequent research relating BIRADS classification and BC 
risk indicating risk due to density and increasing weight should be considered 
independently. Some reports use normalized weight-height relations such as body 
mass index (BMI) or Quetelet"s index (Ql) . In general most studies indicate that 
women with high-risk patterns or increasing densities are more apt to have a lean 
body build. 

[0058] 

Earlier work by Gravelle et al indicate no relation of height with pattern, but when 
stratifying by menopausal status show the pattern-weight relationship is significant 
for postmenopausal women with marginal significance for pre-menopausal women. 
Other work supports the lack of association of pattern-height. Brisson et al indicate 
with increasing weight there is a considerable reduction in high-risk patterns for all 
height categories, but the inverse height association (when controlling for weight) is 
somewhat weaker and not as regular, although visible. Similarly, other research finds 
no dependency on height and pattern. Other work shows that when controlling for 
height, heavier women are more apt to have less dense patterns, and when controlling 
for weight, taller women are more apt to have dense patterns. Kaaks et al found none 
of the anthropometric variables (height, weight, BMI or hip circumference) to be 
significantly associated with BC risk in both pre and post menopausal women in work 
that did not consider tissue characteristics; however, a positive association of BC risk 
associated with waist to hip ratio was observed (RR of 2.79). Similarly, other work 
(without regard to pattern) found no association between BC risk and height. Brisson 
et al shows the percentage of the breast showing nodular or homogenous densities is 
associated with height and weight (positively and negatively respectively); a similar 
trend holds good for homogenous densities with respect to weight but varies little 
with respect to height. Grove et al found women with more dense patterns often have 
lighter body build with or without adjustments for height. The weight-pattern 
association is also supported by a study of autopsy data showing that obesity and 
breast size are good predictors of radiolucent patterns. Longitudinal analysis suggests 
not only are women with low risk patterns more likely to be heavier (large BMI), 
women with high risk patterns that are heavier (large BMI) are more likely to 
experience a favorable pattern change in time. Boyd et al show women with increased 
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radiographic density have significantly higher levels of high-density lipoprotein after 
adjusting for confounding effects, and conclude that differences in fat metabolism 
exist between women with and without increased density. 

Obesity during the pre-menopausal years reduces BC risk due to an associated 
reduction of progesterone. This protective effect diminishes after menopause and 
gradually reverses to increased risk due to bio-available estrogens. 

Family history (FH) is a known risk factor for breast cancer. As stated previously, 
for incorporating risk factors with tissue analysis, serial or non-time dependent risk 
modeling, and CAD nearly uncorrected (or more strongly independent) factors of risk 
bring more power to the decision making process. 

A case control study suggests a similarity in parenchymal patterns in twins, more 
prevalent in monozygotic than dizygotic, implying a familial or genetic influence on 
patterns. This is supported by subsequent genetic segregation analysis that studied 
women with FH, where a significant correlation in sister-sister breast density was 
found when controlling for (or not) many other risk factors, although no significance 
in the mother-daughter densities was found. 

[0062] Research that compared normal women with family history to normal women 

without history shows no difference in patterns. Berkvist et al show, with univariate 
analysis stratified by age or with multivariate analysis, there is no significant relation 
with Wolfe patterns and FH, which is in agreement with subsequent research. Boyd et 
al found a weak association (univariate analysis) with FH and increased density, and 
when considering FH and parity in joint analysis found a weak relation, which agrees 
with subsequent research. Similarly, other researchers found a weak association with 
FH and Wolfe patterns with univariate analysis. In subsequent work, it was shown that 
women with both FH and with greater density proportions are at a significantly 
increased BC risk compared to the elevated risk of FH alone. Saftlas et al also found 
that the combination (greater than the sum) of high risk Wolfe pattern with FH 
exceeds the risk of either effect individually and conclude the two are independent 
influences on risk. This work also shows that women with Nl patterns without FH are 
at minimal (no risk) BC risk. Brisson et al show adjusting for FH (among other risk 
factors) has little effect of the risk estimates associated with parenchymal pattern. 



[0059] 



[0060] 



[0061] 
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Subsequent work shows the means of homogenous and nodular densities are not 
associated with FH in postmenopausal women, and are weakly associated in pre- 
menopausal women. Moreover, other evidence also shows that the risks due to FH and 
mammographic pattern are additive indicating (a) that they are nearly independent, 
and (b) the etiologic role of the two are not closely linked. Women with FH had a 
relative risk of 1 .43 for BC compared to women with no family history of BC while 
women with high risk patterns and FH of are at a relative risk of 4.4 for BC compared 
to women with Nl patterns and no FH of BC. Byrne et al also found the risk due to FH 
and mammographic density or Wolfe patterns are independent which agrees with 
earlier work. 

[0063] Other researchers also found a general lack of association between FH and pattern 
with multivariate methods. However, one report, that investigated the association with 
mammographic density and FH in Japanese women, among other factors with 
multivariate analysis, indicates an association with dense area and FH. Other 
univariate analyses (or similar methods) find little, or no, association with pattern and 
FH, and yet others do find some associations. Wolfe et al also found some relation 
between pattern and FH in a matched case control study. The conclusion drawn from 
this work is that FH and mammographic attributes are practically independent sources 
of BC risk. 

[0064] Understanding the corresponding correlation of breast composition with 

histological features as observed from radiographs are important pieces of the 
pattern-risk premise; this is particularly true when considering that women with 
biopsy-proven benign breast disease show an increased risk for subsequent invasive 
cancer associated with proliferative (epithelial) disorders as compared with non- 
proliferative disease (fibrosis)O 39). Early work by Wellings and Wolfe shows a 
considerable association with histological findings and pattern (a) the low risk pattern 
lacks nodularity, linear densities, or confluent densities, (b) the PI pattern shows 
increased periductal and perilobular fibrosis that is probably related to linear 
densities, where the intralobular connective tissue and heavily collagenized 
extralobular connective tissue are seen in normal lobules, (c) the P2 pattern shows 
increased fibrosis, atypical lobules, and moderate atypical hyperplasia, and (d) the DY 
pattern shows marked fibrosis that is more homogenous than the PI or P2 patterns. 
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To a lesser degree, other work shows some correlation with pattern and histology 
where lobular epithelial hyperplasia and ductal papillary hyperplasia are associated 
with the P2 pattern. This is supported by other work based on autopsy data indicating, 
when stratifying by age for women over 50 years of age (only), there are significant 
associations with high risk patterns and marked intraductal hyperplasia and micro- 
calcification of the terminal duct lobular epithelium. 

[0065] However, earlier research found little (no) association with Wolfe pattern and 
histological findings that indicates the histogenesis of mammary response is not 
related to epithelial alteration in the ductal system. Similarly, other evidence indicates 
no association with histological finding and breast pattern, where the variation in 
Wolfe patterns is attributable to the ratio of fibrous to adipose tissue with no relation 
to the amount of epithelial content. 

[0066] Research based on the percentage of dense tissue in women 40-49 years of age 
shows that women with extreme proportions of density (+ 75 %) are at an elevated 
risk for cancer in situ, atypia hyperplasia, hyperplasia without atypia, and developing 
non-proliferative disease indicating a connection with high risk histological changes 
in breast epithelium. Although, this work indicates a lack of association between 
mammographic density at the biopsy site and histological changes in epithelium that 
might be explained if mammographic densities are caused by changes in the breast 
stroma and not epithelium. This suggests that collagen is the histological feature 
related to breast density detected by radiography that is in agreement with earlier 
work. As expressed in a commentary given by Feig, although dense breast represents 
large amounts of fibroglandular tissue, it is not clear if this tissue is predominantly 
epithelial or fibrotic. 

[0067] As discussed by Madigan et al, roughly 40 of the breast cancers in the U.S. are 

attributable to the known BC risk factors, which was investigated without considering 
density. It would appear plausible to include a tissue related metric in the risk 
analysis. 

[0068] Tne ambiguous nature with age of menarche is in agreement with earlier work in 
developing a comprehensive multivariate model for BC risk prediction that includes 
Wolfe pattern risk. This work indicates that the contribution due to age of menopause 
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is significant but difficult to interpret, age of menarche has a relatively small 
influence, and that menstrual history is better represented as years of menstruation as 
noted earlier. In the comprehensive risk model developed by Whitehead et al Wolfe 
pattern, age of first birth, nulliparity, family history (BC in mother), history of benign 
biopsy, log(weight-lOO), a surrogate for weight, and years of menstruation were all 
found to add significantly to the predictive power of the model. The work showed that 
the interactions with pattern and the other factors in the model were not significant 
indicating each factor is a distinct source of risk. In subsequent work the non- 
interaction of pattern with the other factors was verified; adjusting for differences in 
the non-mammographic factors produced only trivial changes in the BC risk due to 
patterns, indicating pattern is an independent risk factor. Other research supports the 
length of reproductive years as a BC risk when controlling for mammographic density 
in addition to other risk factors such as age at first birth, history of biopsy, education, 
and family history. This work also indicates that (besides weight) adjusting for the 
other factors changes the magnitude in the risk with respect to Wolfe pattern or 
percentage of densities very little, again indicating distinct sources or risk. However, 
earlier work in developing a multivariate BC risk model, where pattern was not 
considered, found that age of menarche, family history, age at first birth, and log 
(etiocholanolone excretion) added to the predictive powers on the incidence model 
while breast feeding, parity and personal history of breast disease did not contribute. 
Likewise other work that investigated the joint effects of reproductive factors and 
body shape suggests the protective effects of late menarche and early first birth are 
greater for pre-menopausal women and both nulliparous and lean women are more 
apt to experience a protective effect due to early menopause. Other work also 
suggests different menopausal status related risks for in situ cancer as related to BMI. 

[0069] The Gail model which includes age at menarche, age at first live birth, number of 
previous biopsies, and number of first degree relatives with BC but not pattern or 
tissue related measures, can also be applied to predict the BC risk over some specified 
interval. More recently, validation research of the Gail model indicates that the method 
can predict the number of cases in a specific risk stratum but has less discriminating 
power on the individual level prediction. 

[0070] The degree of mammographic density is most likely related to the amounts of 
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epithelial cells or tissue in the breast, and if the initiation of cancer is a stochastic 
event, the more cells at risk, the greater the potential for BC incidence, which 
supports the total density area approach, which is supported by the breast reduction 
work. This also suggests that both ipsilateral views should be considered when 
assessing risk. It is very easy to conceive of a geometric figure that is dense in one 
projection and not in the other (for instance a disk). Moreover, maybe it is the density 
acquired from all four views, combined in some fashion, that provides a more accurate 
risk measure. 

[0071] There are degrees of density (due to the extended volume) that are as important 
as the binary measure, and the degree of density (in some crude fashion is related to 
the volume) may be a preferred figure of merit for risk assessment as opposed to the 
dichotomous density labeling. For example, within one image an extremely bright 
area may correspond to densities that extend throughout the breast volume to a 
greater degree (more x-ray blocking ability) as compared with another area that is 
half as bright, but not fat. A texture measure may be deployed in addition to the 
amounts of density in order to enhance the discrimination of the associated risk since 
the P2 (nodular densities) and DY (homogenous densities) breasts appear to behave in 
differently. 

[0072] Tne p resent invention includes methods of modeling mammograms with 

parametric statistical methods as well as non and semi parametric methods. This 
advancement in the art is based on multi-resolution wavelet analysis, spectral 
characterization, and fractal modeling. The present invention includes a risk model 
based on parameters that are not correlated with density measurements as well as 
those that are partially correlated with tissue metrics in conjunction with density 
metrics. The invention (a) incorporates this model into the CAD detection methods, in 
particular for variable threshold adjustments, (b) segregating images out of the 
processing pool that do need to be machine processed which may have a profound 
impact on detection FP rates when the (CAD) algorithms are applied to high risk 
images only, and (c) develops a comprehensive risk model to aid the mammographer's 
decision process that includes measurable image factors such as density as well as 
other risk factors. The evidence indicates that a starting place for this analysis will 
include FH, length of reproductive years, and a binary variable for parity, in addition to 
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density measurements. The modeling will include making assessments of risk today 
for a particular subject without regard to the past as well as considering compounding 
risk over time as in the Gail approach: that is, assessing the absolute risk of breast 
cancer over some specified time interval. 

[0073] Aging in itself is a risk factor. This follows from noting that the breast cancer 

incidence curve is an increasing function of age. Even though the incidence risk is an 
increasing function of age, the rate of change slows (still positive) after the age of 50 
indicating that breast cancer is hormonal in nature, since 50 is roughly the average 
age of menopause, which is discussed in more detail below. The inference from this is 
that breast density (or pattern) will shift to a lesser degree (tracking risk) with age as 
Wolfe noted early on. 

[0074] Here the longitudinal tissue pattern change is discussed. Among the various 
reports there is almost uniform consensus that the pattern distributions have a 
tendency to regress from high risk to low risk with age. Additional evidence, obtained 
from autopsy data of non-selected non-referred women between the ages of 1 5 and 
98, indicates that age and postmenopausal status are very strong predictors of 
radiolucency. However, one report shows little difference in density (average 
measurements) in a cohort stratified at 50 years of age, although with considerable 
variation. 

[0075] There is considerable variation in this analysis, and the interpretation of the data 
must be taken carefully. For instance, cross-sectional (or case control) studies 
stratified by age give a static view of the pattern distribution. The inference drawn 
from the cross-sectional work is that the time evolution of the breast tissue is the 
mechanism driving the age related distribution. A stronger method for determining 
temporal pattern change is derived from the serial studies that track a group of 
women over extended time intervals with initial and final mammography assessments. 

[0076] in general, the cross-sectional work demonstrates that there is a considerable 
tendency for the more dense patterns to regress to less dense patterns with the DY 
proportion showing significant shifts with increasing age with regression occurring 
within the DY group. Other work, to varying degrees, shows similar trends in 
decreasing dense patterns with age. The average age for women with P2DY patterns is 
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about 53 years, and about 58 for the low risk Nl PI group. Stomper et al show a 
progressive decrease in percentage density with age, where 39% of the women in the 
25-29 age group have 90% or greater density, whereas 76% of the women in the 75- 
79 age group have predominately fatty type patterns with only 6% of the women 
having 90% or greater density. 

[0077] When stratifying the age analysis by type of density and menopausal status it 
appears that homogenous densities and nodular densities behave differently. The 
mean percentage of the breast showing nodular densities decreases with age 
regardless of menopausal status whereas the mean percentage of the breast showing 
homogenous densities decreases sharply with age for pre-menopausal women only 
and then stabilizes. In subsequent work Brisson et al verify the dense tissue reduction 
with age and show that for all age groups postmenopausal women have less dense 
tissue. Likewise, multivariate work shows that age has an opposite effect on pre- and 
post menopausal women in that pre-menopausal women age 40 and above have twice 
the risk compared with women under age 36 while post menopausal women over age 
55 have half the risk of postmenopausal women aged under 51 years of having P2DY 
patterns. This is supported by research suggesting the P2 pattern is more common in 
older (or post menopausal) women or younger (pre-menopausal) women. Inferences 
drawn from earlier work by Grove et al implied menopausal status was more 
important for predicting pattern than age. However subsequent work weakly 
confirmed this association. 

[0078] only a few studies have examined the serial pattern change with longitudinal 
analysis. One report, which tracked the subjects over a mean interval of 1 5 years, 
indicates important points (a) the DY pattern regresses to Nl , PI , and P2 (in 
ascending frequency) at mean ages of 55, 57, and 42, respectively, and (b) small 
percentages of the P2 patterns regress to PI (only), with the proportions of P2 
increasing with time, and (c) the PI patterns show no change and (d) a considerable 
proportion of the Nl breasts progress to PI such that the Nl patterns show a small 
decrease over time. Similarly, other work shows (1 2 year interval) the largest 
regression from dense parenchyma occurs in cohorts aged 45-54 years and 35-44 
years of age. Although the research shows that the proportions of P2 patterns are 
somewhat stable over many age groups, there is a considerable change in patterns 
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that were originally designated as P2 to less dense patterns, and the PI shows 
dynamic behavior in both directions over all age groups. This is consistent with earlier 
work indicating the P2 patterns are relatively stable in time (7 year interval), although 
some do regress to less dense patterns. This work also divides the DY group into 
three subgroups of increasing density and indicates considerable regression within 
the DY group associated with the most dense DY patterns; about one third of the 
P2DY patterns regress to a less dense state. Salimen et al also show a considerable 
number of P2 and DY patterns changing in time (8 year follow-up) with women in the 
45-53 age group almost twice as likely to regress as compared with women 40-44 
years of age. 

[0079] For all practical purposes, the reports indicate that the high-risk patterns regress 
with age. Since the commencement age for screening is often between 40 and 50 
years of age it is difficult to describe this dynamic process with great precision. It 
should be understood that the evidence presented here indicates no clear lines of 
demarcation such as age of menopause or some particular chronological age where 
the patterns change abrubtly; any apparent boundaries are an artifact of the grouping 
or study design. This work indicates that any attempt to successfully apply automated 
change detection (serial analysis) to digital mammograms will have to cope with time- 
changing densities, which for the most part are normal occurrences. This is in 
contrast to the often-stated problem that differing x-ray projections associated with 
the temporal views are an obstacle for serial analysis that must be accounted for prior 
to analysis. 

[0080] involution is the aging process of the breast, where by the functioning elements 
are replaced by fat. There are relatively few studies that address this mechanism 
directly, even though the results presented in the previous section are due to 
involution effects. This process, which begins many years prior to menopause, can be 
characterized in two stages (a) pre-climactic characterized by moderate atrophy of the 
glandular epithelium and a decrease in the amount of lobular and acinar tissue, and 
(b) the menopausal stage, characterized by a progressive loss of all glandular tissue 
accompanied by proportional increase of fat and connective tissue. It is thought that 
this process is hormone-related. Related work, where mammary glands were collected 
at necropsy from non referred women, shows (a) the proportion of epithelial tissue in 
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the breast is a maximum in the third decade of life (pre-menopausal), which 
decreases afterwards, (b) the average lobule volume decreases with age, while the 
lobule volume proportion of epithelial tissues increases until around 50 years of age, 
and (c) the fibrous tissue volume proportion declines until about 50 with the trends in 
these tissue types reversing after this. 

[0081] The breast cancer incidence vs. age curve (age specific incidence rates) increases 
with age, but has an inflection point at about 50 years of age, where the increase 
slows (but still is positive) appreciably with older age. This curve is consistent with the 
hypothesis that the loss of glandular tissue leads to a reduction of risk; that is the risk 
is still increasing but the increase has a slower rate of change. This plot consists of 
two linear segments with differing slopes that interact at around 50 years of age 
(about the average age of menopause). For most (non hormonal related) cancers there 
is a linear log-log relation in the incidence curve (without the kink). This, taken in 
conjunction with certain BC risk factors, may indicate that breast tissue does not age 
in the conventional manner instep with the calendar. 

[0082] The rate of involution may be the true variable associated with risk. The 

associations of BC with hormonal interactions and involution are suggested by the 
connection with high and low risk patterns. Oophorectomy appears to have a 
protective effect against BC that may be due to artificially induced involution, whereas 
early menarche seems to contribute to enhanced risk. Likewise, late pregnancy, which 
is associated with the prevalence ofdense patterns, can be considered delaying the 
involution process; thus increasing the risk for BC due to the attendant increased 
proliferation in the ductal-lobular alveolar system. This suggests that factors that 
inhibit involution or induce an extended life span of the functional elements m the 
breast at older age may contribute to increased risk. However, there is evidence that 
indicates the amount of epithelial tissue in the breast and number of pregnancies is 
not related. Similarly, a discussion related to aging and breast development provided 
by Russo et al indicates that the initiation of cancer requires the interaction of a 
carcinogen with an undifferentiated and highly proliferating mammary epithelial. 
These authors also discuss pregnancy and how the post partum breast contains more 
glandular tissue than if the pregnancy never occurred. 
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[0083] The concept of early involution providing some protective affect from BC is 

supported by work showing that Native Americans shift to less dense patterns at an 
earlier age compared with Caucasian women, and are also at reduced risk of BC 
compared with Caucasian women. Likewise, as discussed below, other influences such 
as HRT may delay involution, and thus induce an additional risk. 

[0084] Involution may not necessarily occur uniformly across the breast. The related 

histopathological description is provided in detail elsewhere. Roebuck et al provides 
another pattern classification, which parallels Wolfe"s work, that considers dysplasia 
(density) and prominent duct opacities as shadowing artifacts on the normal breast 
structure. This classification accounts for severity, extent, activity, and accommodates 
ageing. This study discusses how the involution process affects the classification and 
appearance of the breast, and that P (related to ductal components) or D (related to 
dysplasia) shadowing can only be considered normal in areas where glandular tissue 
remains. Roebuck et al also state that once P shadowing has become established, no 
reduction with time occurs, although it may become more pronounced with age. 

[0085] Russo et al discuss lobular development (three types of lobs), and transitions from 
one to the next. This work indicates (a) nulliparous breast parenchyma is made up of 
undifferentiated structures such as terminal ducts and lob 1 , which remain the 
predominant structure, whereas (b) history of full term pregnancy correlates with an 
increase of lob 3 (the proportions of lob 3 are significantly lower since they are 
occasionally seen). These researchers also indicate that (a) lob 2 is present in early 
years but sharply decreases after 23 years of age, (b) menopause shows an increase of 
lob 1 and associated decreases in lob 2 and lob 3, (c) at the end of the fifth decade 
the breasts of nulliparous and parous women are primarily composed of lob 1 , and 
put forth an interesting premise: the endocrinological and physiological changes 
triggered by pregnancy in prior years might have imprinted permanent biological and 
or genomic changes in the breast that may impact the BC development potential even 
though they might not be manifested at a morphological level. 

[0086] 

Ductal carcinoma is the most frequent breast malignancy and has been identified 
as originating in the terminal ductal-lobular unit (or lob 1). Likewise, lob 2 has been 
associated with this malignancy development, but no association with lob 3 and the 
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development of malignancies has been found. Experimentally induced mammary 
cancer research on tissue derived from reduction mammoplasty surgeries shows in 
vitro cells derived from lob 1 and lob 2 express more readily, changes that are 
indicative of neoplastic transformation suggesting that these cells are more 
susceptible to the transforming effects of genotoxic agents. 

[0087] These ideas suggest that densities measured in the nulliparous women are not 
equivalent to densities measured in parous women. This is consistent with work 
showing the loss of dense patterns due to parity primarily influences the P2 
population and that the P2 densities have an extensive nodular component. 

[0088] The masking effect is an argument used to explain (or better refute) the pattern- 
risk premise when assessing incidence risk. The thought behind masking is that 
developing abnormalities are sometimes missed or obscured in dense fibroglandular 
tissue, resulting in delayed detection. This should (a) decrease the prevalence rates 
with respect to dense patterns and (b) artificially inflate the incidence risk estimates 
for dense patterns in the short term because abnormalities missed will manifest 
shortly thereafter, which would then cause the risk to fall over extended intervals. If 
masking is in effect, evidence that shows decreased screening sensitivity due to dense 
breasts, in part, provides corroborating evidence. Similarly, examining the factors 
associated with missed cancers, which are closely related to sensitivity, should help 
shed light on this subject. However, there may be confounding effects such as age 
dependent tumor growth. Rapid growth and missed cancers are closely related 
through a category of abnormalities termed interval cancers. 

[0089] There is a considerable amount of research devoted to understanding the 

characteristics of both missed cancers and closely related interval cancers. Missed 
cancers are related to sensitivity in that they are major factors associated with reduced 
screening sensitivity. In some instances missed cancers can be considered as a subset 
of interval cancers. That is if a cancer is currently diagnosed in between screening 
visits and observable in retrospect in the previous mammogram, it is considered as a 
missed cancer due to an error in reading. A missed cancer is also a definition used 
when a cancer is detected in a current screening mammogram that is also visible in 
the previous screening image with retrospective analysis. 
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[0090] Interval cancers are cancers that are detected, in some fashion, in between regular 
screening visits. If the previous screening mammogram, reviewed in retrospect, is 
negative, the cancer can be considered as either a true interval surfacing cancer or 
mammographically occult at time of initial screening. In this sense, mammographically 
occult implies that the tumor was present but under the sensitivity limits of 
mammography. 

[0091] Some authors explicitly discuss a class of cancers termed true interval cancers, 
and yet others regard this same situation as occult cancers. True interval surfacing 
cancers and occult cancers following the definition above are indistinguishable, and 
often lumped together. The characteristics of these cancers (true interval surfacing 
cancers) are important in terms of masking. For instance in the limiting case, if all 
interval cancers are occult cancers that were present previously but not 
mammographically observable, and some known proportion of occult cancers occur in 
dense breasts, the odds ratios in risk analysis can be adjusted to find the limiting risk 
ratio relative to the less dense patterns. These definitions, as associated with interval 
cancers also apply in an analogous way for screen detected cases with retrospective 
image analysis. 

[0092] Another definition of occult follows from a cancer that is currently detected by 

means other than with mammography that is not visible in the current mammogram, 
such as palpable lumps. If these cases show some breast pattern dependency, a 
circumstantial case can be made for masking, since it may be safe to assume that 
these cancers were present before becoming symptomatic and were not visible. These 
cases are referred to as occult at diagnosis and truly define the lower limits of 
mammographic sensitivity. 

[0093] A closely related case is when subtle signs are found in retrospect, which 

generates another classification (sort of a hybrid interval classification of occult-error 
and interval developing). Some studies discussed below investigate interval cancers 
specifically, and other studies investigate previous mammograms of cancers detected 
presently at regular screening using the same classification as discussed above. 

[0094] The masking effect is partially supported by the incidence studies that show the 
long term risk diminishing for the more dense patterns following initially elevated 
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risks. However, there may be other reasons that provide equally valid explanations for 
the rising and falling risk estimates as discussed below, van Gils et al studied the 
masking effect and show elevated incidence risk for women with dense breasts (25% 
density or greater) up to the 3-4 year interval and then decreasing risk in the 5-6 year 
interval, which could be interpreted as masking to some degree; it is important to 
note that the density determination was made at the initial screening. Work by Sala et 
al also show rising and falling incidence rates, where it is suggested that masking and 
accelerated cancer growth rate influences may be operating in tandem. Whitehead et 
al provide an argument based on a hypothetical screening program, where it is 
assumed equal incidence probabilities for P2DY and N1 PI patterns that indicates the 
masking effect has a minimal influence on inflating the BC risk and does not explain, 
completely, the elevated risks associated with high risk patterns. This can be 
understood by considering that masking is always present, and hence the delayed 
detection in the screening environment is not affected because there is a continuous 
folding back of cases due to delayed detection. 

[0095] 

Sala et al show the risk for interval cancer is higher for women with P2DY patterns, 
in particular the DY group, compared with screen-detected cancers. Subsequent work 
indicates that occult cancers occur more frequently in dense breasts and younger 
women with the two posing a greater risk than either alone, suggesting independent 
influences. When comparing the subset of invasive ductal cancers, women with P2 and 
DY patterns are at a significantly elevated risk of having a grade three cancer 
compared with women with less dense patterns, which may be interpreted as rapid 
cancer growth. Subsequent work (case control study) by Sala et al indicates the odds 
of having invasive cancer compared with cancer in situ are higher for women with 
P2DY patterns, and that these women are more likely to have lymph node 
involvement. Here again the interpretation is either masking or possibly a dynamic 
process associated with dense breasts that could be responsible. Similarly, other work 
shows an increased risk for women with P2 similar patterns is more pronounced for 
women with poorly differentiated more aggressive grade three cancers. Work by Adler 
et al also indicates that the number of false negative cases (FN) is almost twice for 
invasive lobular carcinoma, which is cancer of epithelial in origin compared with 
ductal carcinoma with FNs found more often in glandular followed by fibrous and then 
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by dense breasts and in younger age groups. 

[0096] Evidence indicates decreased sensitivity in mammography screening associated 
with dense patterns, although not unequivocal, which again renders support for the 
masking hypothesis, although risk assessment is not addressed in this work. Methods 
used to calculate mammographic sensitivity are discussed elsewhere and in references 
therein. Work by Bird et al indicates that missed cancers, which include cancers 
incorrectly diagnosed at reading and cancers visible in previous mammograms, are 
more likely to occur in dense breasts, in addition to showing cancers that present as 
developing opacities are also a cause for missed reading. Similarly, other work also 
shows missed cancers (error in reading) are significantly lower in mammographic 
density, relative to the surrounding tissue, and often only seen in one of the two 
views. Kerlikowske et al show that the sensitivity of mammographic screening is 
highest for women over 50 years of age who have primarily fatty breast density, 
whereas for women under 50 years of age breast density did not influence sensitivity. 
Similarly, evidence also shows age related sensitivity measures without regard to 
pattern, where sensitivity increases with age. Rosenberg et al show that age per se for 
women over 40 is a minor determinant of sensitivity, but indicates breast density is an 
important contributor for decreased sensitivity in particular for women 40-50 years of 
age contrasted with older women. Other work indicates that the poorer sensitivity in 
younger women may be due to rapid tumor growth, and shows no decreased 
sensitivity for women 40-50 years of age. Work by Lehman et al indicates that women 
with extremely dense breasts are almost twice as likely to have a false positive 
mammogram independent of age. Ma et al show that breast density, smaller tumor 
size, and lobular histology are independent factors that contribute to the false 
negative mammogram and indicate that biological factors are the root cause for 
misreading. Evidence also indicates age-related differences in patients with non- 
palpable breast carcinomas in a biopsy study of occult tumors, where the majority of 
tumors were invasive ductal carcinomas with median patient age of 59 years. Women 
with ductal carcinoma in situ had median age of 50 years whereas women with 
invasive ductal carcinoma without associated intraductal tumors had median age of 65 
years. 

[0097] Without regard to age or mammographic signs, evidence indicates missed 
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cancers, including screen and interval detected cancers, account for roughly about 1 2 
% -30 % of the cases. These figures of merit are supported by earlier work that the 
shows roughly 30% of the missed cancers are due to observer error, and about 1 /3 of 
the missed cases show some subtle signs in retrospect. Yet other work indicates that 
roughly 57 % of the interval cancers are truly interval surfacing. Porter et al indicates 
that the true interval cancer rate is about 1 7%. Moreover, there is evidence indicating 
that interval cancers (all types) are found more often in younger women compared 
with screen detected cases. A review of some earlier work in this area can be found 
elsewhere. 

[0098] With respect to pattern, research comparing missed cancers, occult cancers 

(interval surfacing), and occult at diagnosis cancers detection indicates that there is 
little association between the density groups (PI Nl and P2DY) or type of detection 
category when combining screen detected with interval detected cases. Burrell et al 
using a slightly different detection classification and studying strictly interval cancers 
indicate a similar lack of association; however, all the occult at diagnosis cases were 
found in the P2DY group. This is in agreement with earlier work showing occult at 
diagnosis cancers were more likely found in women with low Ql, and with women with 
P2DY patterns compared to women with screen- detected or interval detected 
cancers. This work indicates that about 1 6% of the interval cases are occult at 
diagnosis (but included cases with subtle signs), and finds no pattern dependence for 
missed interval cancers. Similar associations with high risk Wolfe-patterns and occult 
cases of this nature were also found in earlier work. One report that investigated 
interval cancers with respect to BI-RADS tissue composition finds women with dense 
breasts are at a significantly high risk for interval cancer when controlling for other 
factors compared with screen detected cases, but showed that when stratifying 
interval cases by age, women between 40-49 years of age with dense breasts were 
not at an elevated interval risk when compared to women with less dense breasts. 

[0099] Tnjs work j nc ji cates that masking due to dense breasts is a factor to consider, but 
the magnitude is not that great. However, this should be qualified with the 
understanding that much of the work investigating interval, occult, or missed cancers 
does not stratify the analysis by pattern and age simultaneously. In addition, the lack 
of analysis using proportions of density (rather than Wolfe patterns) results in less 
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than a complete picture. One safe assumption is that a well defined incidence-risk 
study should make adjustments for cases that are missed in the previous 
mammogram. 

[0100] The variation in the sensitivity-age relation may be in part due to exogenous 

hormone use. Evidence indicates that estrogen replacement therapy, in conjunction 
with increased breast density, are significant factors that contribute to decreased 
sensitivity in women over 50. Similarly, other research indicates reduced sensitivity 
and specificity (false positive or FP rate) for women over 50 years of age using 
estrogen replacement. However one report that studied the effect of HRT in current, 
previous, or never users indicates no decrease in sensitivity due to HRT with a 
marginal decrease in specificity. However, other work indicates that women on HRT 
are almost at twice the risk of experiencing an interval cancer relative to non-users, 
which is interpreted as a sensitivity decrease. 

[0101] Cancer dynamics and tumor growth rates are important factors to consider when 
designing screening systems. These factors may also be related to BC risk and 
reproductive factors and are discussed below. The time that a particular disease 
remains in the pre-clinical stage is termed the sojourn time, and is an important time 
interval to consider when assessing the efficacy of a screening program. Often this is 
replaced by the mean sojourn time (MST). Likewise, an important related parameter 
for assessing tumor growth rates is the tumor volume doubling time (DT). 

[01 02] The sojourn time is related to the screening sensitivity and has an influence on 

screening practice. Assume that a disease is detectable by some imaging modality at 
time TO, but is asymptomatic, and that the disease exhibits symptoms (becomes 
clinically apparent) at time Tl , the interval Tl -TO is the sojourn time. This is the total 
time that would have elapsed when the disease could be detected without diagnostic 
intervention. Now assume that the individual is screened at time T2 (TO < T2 < Tl) 
with the correct pre-clinical stage diagnoses, the lead time or the time that the 
disease is brought forward is Tl -T2; the probability that the diagnosis is correct in 
this interval is the sensitivity of the test. The sojourn time provides an absolute upper 
limit to the lead-time obtainable and gives an indication of the effective screening 
interval. The shorter the sojourn time the shorter the screening interval. 
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[0103] There is considerable evidence indicating tumors progress faster in younger 

women from the preclinical to clinical phase, implying faster growth rates and shorter 
sojourn times compared with older women. Other related work indicates similar 
effects. Using a Markov-chain model to estimate tumor progression rates, research 
shows shorter sojourn times for women between 40-49 years of age compared with 
older groups of 50-59, and 60-69 years of age with sojourn estimates of 2.44, 3.70, 
and 4.1 7 years, respectively, which is consistent with earlier work. This is also in 
agreement with earlier work showing similar trends in age dependency lead times. 
Although, Brekelmans et al provide even shorter age specific sojourn times of one 
year for the 40-49 age group. 

[0104] 

Related research shows that the proportion of tumors capable of dedifferentiation, 
that is the ability of a tumor to change grade, appears greater for the 40-49 year age 
group. This is in agreement with work discussed earlier suggesting true interval 
cancers (aggressive pathologies) are more frequent in younger women, are larger in 
size, and may be biologically different compared with cancers detected at regular 
screening. Earlier research also indicates that interval detected cancers may be 
biologically different from cancers detected at regular screening; that is they are fast 
growing, and more often found in younger women, which is in agreement with other 
research. DeCroote et al also suggest interval cancers are a more aggressive form of 
cancer compared with screen-detected disease. Gilliland et al show that (a) 
mammography is more likely not to detect cancers with high proportions of 
proliferating cells (Ki-67), (b) interval cancers have a greater potential for dis- 
regulation of cell cycle (a potential for genetic instability) as measured by p53 
expression when comparing interval with screen-detected cancers, and (c) younger 
women have higher proportion of rapidly proliferating and aggressive tumors. 
Similarly, other evidence indicates that the mean sojourn time for women in the 50-69 
age group is dependent on the tumor histology, ranging from 2-7 years, whereas 
women in 40-49 age group show no such histological dependency with a mean time 
of about two years. Subsequently, Tabar et el discuss the evidence that indicates 
tumor grade, size, and lymph node involvement progress more rapidly in younger 
women, which is in agreement with similar research. In contrast, other work also 
indicates that interval cancers are more likely to emerge in women between 40-49, 
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but shows little evidence of biological aggressiveness. 

[01 05] Other work shows that the tumor volume doubling time is appreciably faster in 

younger women. It should be clear that both sojourn times and tumor doubling times 
are slightly different measures, but are closely related approaches. Work, using an 
exponential growth model approximation, indicates age dependent volume DTs of 
about 1 57 days for women between 50-70 years of age and 80 days for women 
younger than 50. The age related doubling times are similar to those estimated 
earlier. It is interesting to note that applying the exponential tumor growth model, 
and assuming malignant growth starting from a single cell, or small cell clusters, 
suggests most breast cancers start growing around the age of 40. One report 
indicates no relation between tumor growth rates and Wolfe patterns. 

[01 06] There are three interacting factors that are responsible for the confounding 

pattern-risk incidence studies (a) younger women are more apt to have dense breasts, 
(b) evidence indicates that malignancies progress faster in younger women, and (c) 
mammography is less sensitive to the more dense breasts. There is some research 
indicating that there may be an association with pattern and cancer dynamics. Ciato et 
al found an increased prevalence of stage one tumors in the less dense images which 
is supported by other work showing women with P2 or DY patterns are more likely to 
have grade three cancer as opposed to women with Nl or PI patterns that are more 
likely to have tumors characterized as grade one or grade two. 

[0107] 

The term parenchymal density refers to a composite of stromal and epithelial 
tissue. The evidence indicates that there is some interplay between these constituents, 
and the stromal content may act as a growth factor reservoir that not only helps 
promote tissue growth but also plays some regulatory role in cancer growth. The 
stromal- epithelial growth factor associations and the sex steroidal relation with 
mammary tissue proliferation are provided elsewhere. Since the evidence indicates 
that cancer progresses faster in younger women, and younger women are more likely 
to have dense tissue, it would be reasonable to make the inference that there should 
be a density dependent tumor growth rate exhibited, particularly in younger women, 
which would then provide support for the pattern-risk premise. One report in this 
area, implemented some time ago, with Wolfe pattern analysis that found no 
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association with growth rates and pattern, which is suggested in another report that 
shows true interval surfacing cancers are not favored in the P2DY combined group. 

[01 08] The conclusion drawn from this survey is that many interactions are at play. If the 
pattern-risk premise is true, then more tumors will manifest in dense breasts. These 
tumors are more likely to grow faster and thus influencing the incidence rates in the 
short term since dense breasts are more common in younger women. It is equally safe 
to assume, that some abnormalities are masked, most often in dense breasts, which 
will emerge and artificially inflate the incidence rates, but to a small degree. 

[01 09] A review of the work in HRT and risk given by Jacobs et al indicate HRT use 
presents an additional BC risk. Moreover, other research in Wolfe pattern risk 
assessment and HRT use shows that only women with DY patterns on HRT are at a 
significantly elevated risk, indicating a synergistic risk affect for the combination. The 
idea here is to understand how this treatment influences the temporal image 
appearance. Here again there are cross sectional studies as well as incidence studies, 
where the pattern is assessed prior to therapy and during therapy with the latter 
giving a dynamic view of tissue change. 

[01 10] In general the cross sectional studies indicate women undergoing HRT treatment 
are more likely to have more dense patterns with one exception. Likewise one report 
indicates increased density for women over 5 5 years of age with no difference for 
younger age women. 

[0111] Many of the dynamic (serial change) studies indicate that there is an increase in 
percent density (proportions of dense tissue) influenced by HRT treatment, with two 
exceptions that indicate the treatment is not responsible for the increased prevalence 
of high risk patterns in women on treatment with no shift toward dense Wolfe patterns 
with the use in cancer patients. Other research indicates that HRT use influences 
unfavorable change in Wolfe patterns with no impact on favorable change. 

[0112] 

The evidence indicates the increasing density shift takes place in the short term of 
approximately 1 -2 years, and the degree of change is dependent on the treatment 
regimen, with the greater influence due to the combination of estrogen and 
progesterone rather than estrogen alone. Similar effects were also found with 
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ultrasound analysis, and magnetic resonance imaging work. 

[01 1 3] Work by Leung et al shows that the change takes place in the short term and that 
no significant density change was noted at the 5 year mark compared with 1 8 months. 
Work by Kaufman et al suggest that women on HRT for five years or longer experience 
a delay in the involution process compared to women not on the treatment that show 
the shift to less dense patterns with age. This idea is supported by subsequent 
research that shows a significant shift to less dense patterns for non-HRT users when 
compared with HRT users, where 80% of the users show no shift by 5 5 years of age. 
Even though the work by Stern et al shows that a small portion of the women 
experience interval density increases with treatment, the work suggests that there are 
significant density decreases in non-users of 55-64 years of age while the density is 
maintained in users, thus supporting the delayed involution premise. Likewise, a serial 
DM study also suggests that HRT use delays involution rather than increasing the 
actual density proportions. The delayed involution observation is consistent with the 
comparisons of I^RT and menstrual cycle breast epithelium morphology discussed 
above. One study, that investigated BIRADS composition classification indicates that 
HRT initiators are more likely to increase in density, discontinuers are more likely to 
decrease in density and continuers are more likely to have increased density, all 
relative to nonusers with older initiators having a greater risk of increasing density. 

[01 1 4] Quantitative assessments of the shift indicate that 73% of the cases experience a 
shift to more density with a mean increase of a little less than 7% with greater shifts 
for those with less dense baseline measurements; this work also shows a shift in 
Wolfe patterns from low to high risk in 24% of the cases. Earlier work by McNicholas 
shows that 27% of the cases showed 1 0% or greater increases in glandular density. In 
addition, two specific cases of emerging densities in the short term due to HRT are 
discussed and shown in detail by Doyle and McLean. More recently, other work shows 
that only 8% of the women experience a density shift. 

[0115] 

Cyrlak et al discuss some of the earlier work in HRT and breast density increase, 
where a pictorial essay of the influence is provided. The authors discuss the 
interaction of progesterone with estrogen and how the two act in synergism in the 
distal portion of the duct favoring differentiation into acini, and promoting lobular 
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growth; and that in diametric opposition, progesterone converts the proliferative 
effects of estrogen on the ductal cells into cellular differentiation. A cross sectional 
biopsy analysis of cell proliferation shows that women treated with estrogen and 
progesterone are associated with greater breast epithelial proliferation and cell 
density than those treated with estrogen alone, or those with no treatment with the 
effects localized to the duct-lobular unit. These researchers also indicate that the 
breast epithelium morphology in the combination treatment is consistent with that 
observed in mammary tissue during the luteal phase in cycling women and is similar 
in amount. In related mouse mammary studies evidence indicates the HRT use 
increases mammary cell proliferation in surgically induced menopausal mice. 

[01 1 6] As with much of the research in this area, no sure conclusion can be made as to 
the real influence of HRT. The compensating variables at play are difficult to isolate. 
The natural aging (decreasing time related density) is most likely different for each 
woman and the HRT involution retardation and or cell proliferation may be 
simultaneous constituents for the HRT-density relation. 

It is quite reasonable to expect that menstrual timing (MT) would have some 
impact on the breast tissue, since natural cycles result from hormonal influences. In 
fact, there is evidence showing MT has an impact on the breast volume size and image 
appearance. 

Magnetic resonance (MR) imaging research shows that there are significant 
differences in total volume, parenchymal volume, parenchymal water content, and T2 
relaxation, but not fat volume or Tl , depending on MT, with the parenchymal tissue 
volume highest in the second half of the cycle, termed the luteal phase. This is 
consistent with subsequent MR work showing increased water and fibroglandular 
content in the later menstrual phase. 

[01 1 9] For reference, the relations with the MR signal and breast tissue in terms of 
radiolucency are: high water (blood), moderate water (epithelial), low water 
(connective tissue), and no water (fat). Similarly, work by White et al shows that a 
smaller proportion of women have extremely dense breasts in the follicular phase 
(first phase) of their menstrual cycle with the effect more exaggerated for women with 
a more lean body build. 



[0117] 



[0118] 
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[01 20] Evidence suggests that dense mammograms are more difficult to analyse, which is 
in accordance with work showing significantly increased chance of a FN mammogram 
assessment when the image is acquired during the luteal phase. A discussion given by 
Baines suggests that menstrual cycle and the related endogenous hormones may hold 
the secret of why younger women do not benefit from screening and subsequent 
therapy. 

[01 21] Menstrual timing is a consideration when acquiring images for automated serial 
analysis. If possible, the images should be acquired within the same phase and 
preferably in the follicular phase for pre-menopausal women. In addition, Spatt et al 
make the observation that if radiation therapy is to induce a carcinogenic mutation it 
would have to do so when the ductal cells are in the G2 phase, and that there are 
greater proportions of cells in this condition in the progestational phase (luteal) of the 
menstrual cycle. 

[01 22] As with other aspects of this analysis there is some evidence that provides counter 
intuitive view. Research, using mammography in conjunction with urine and blood 
analysis to investigate endogenous hormone levels during the luteal phase, shows 
higher levels of estrogen and prolactin and lower levels of progesterone in women 
with Nl and PI patterns as compared with the high risk P2 and DY group. The authors 
suggest that the reverse would be true if endogenous hormones were responsible 
(causal) for dense tissue. 

[0123] There are ongoing debates as to whether or not OC use presents as an additional 
BC risk factor, with many studies finding no association. Recent evidence suggests 
long-term use is coupled with increased risk for older women. 

[01 24] Earlier univariate analysis shows a high proportion of low risk patterns in women 
who have used OC (pre-menopausal). In support, subsequent univariate analysis 
stratified by menopausal status, shows an increase proportion of Nl patterns and a 
decrease of P2 patterns for pre-menopausal women who have used OC in, although 
not significant. In is subsequent multivariate analysis this variable was not found 
important. However, one report indicates no significant association with pattern and 
OC use in the general population, but when stratifying by menopausal status shows 
that postmenopausal women who have ever used OC were more likely to have high 
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risk patterns, de Stavola et al, with univariate analysis, show a positive association 
with OC use and low risk patterns in pre-menopausal women and an inverse 
association for postmenopausal, but with multivariate analysis the associations were 
not found significant. 

[01 25] A cursory discussion of the interplay of hormones, age, breast tissue, and general 
reproductive characteristics might help shed some light on variation in the pattern- 
risk analysis. This may also provide some insight into developing both well designed 
studies and experimental approaches to help resolve some of ambiguities involved 
with the tissue-risk phenomenon. Normal breast development depends on age, breast 
tissue composition, and past/present hormonal environment. As expressed by Pike et 
al and borne out in the work by Whitehead et al, the breast age, either defined as 
starting at menarche, taking in account the distance between menarche and 
menopause, or the distance between menarche and age at first birth are important 
time-related factors to consider, as opposed to the chronological breast age as 
discussed earlier. Boyd et al also provide a similar discussion, premised on this earlier 
work, indicating that the breast tissue, in some way, may provide the cumulative index 
of the past and present hormonal interactions. Work by Rockhill et al also provides 
supporting evidence that the various ways of representing age may interfere with the 
risk analysis. The point to stress here is that it may be worth exploring methods of 
normalizing the analysis to account for age of menarche, and perhaps stratifying the 
analysis according to the distance between menarche and menopause. A good 
discussion of reproductive history, cell proliferation and breast cancer risk is also 
provided by Pike et al. 

[01 26] Much of the density related work is an effort to correlate known risk factors with 
corresponding changes in density (or pattern). If the relationships hold, a mechanism 
then exists to measure if a particular action influences risk; that is breast tissue in 
some way may be considered as a biomarker for environmental or hormonal 
influences. There is evidence indicating that this holds for some important cases not 
discussed above. A case control study shows significant reduction in the incidence of 
mammographic density in postmenopausal cancer patients on tamoxifen therapy (a 
treatment that increases disease free survival rates and mortality for cancer patients). 
As expected the incidence of high risk patterns with the cancer patients relative to the 
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normal control group shows a significant excess of high risk patterns at the initial 
screening. However, the high risk patterns of the cancer group regress in time due to 
therapeutic interaction; the study shows that most change due to treatment occurs 
during the 10-25 month interval, which is in agreement with earlier work. Similarly, 
subsequent reports show a density decrease in high risk women on tamoxifen therapy 
over a two year interval, which is consistent with work by Brisson et al indicating a 
shift to lesser densities (apparent 1 -3 years) primarily for women under 50 years of 
age. These reports coincide with histological studies of pre-menopausal cancer 
patients on tamoxifen treatment showing reduced proliferative activity of mammary 
tissue. Likewise in a study of hormonal contraception that uses gonadotropin agonist 
in addition to low doses of estrogen and progestron, which attempts to minimize 
exposing the breast epithelial to these steroids, shows substantial reductions in 
breast density over a one year follow-up relative to the baseline reading. 

[01 27] In pre-menopausal non-pregnant women, epithelial cells are in a dynamic state 
undergoing a continuing life and death process due to cyclic ovarian (estrogen - 
progesterone interaction) function. In these women, the maximum breast epithelial 
proliferation rate is in the second (luteal) phase of the menstrual cycle, which is also 
suggested by the MR work discussed previously. This repetitive process increases 
probability of somatic cell mutation and thus presents a BC risk. This is consistent 
with the differences in BC incidence curve rate of change in menarche-menopause 
years compared with the postmenopausal period. It is also interesting to note that the 
breast cancer incidence rate of change for women who have undergone oophorectomy 
(artificially induced menopause) below the age of 3 5 years of age parallels that of 
postmenopausal women. Although it is important to note that, to the contrary, earlier 
work shows increased endogenous hormones in the luteal phase, but finds no relation 
with high risk Wolfe patterns and endogenous hormones. 

[01 28] other research indicates that breast cell proliferation is similar over OC induced 
cycles a as with normal menstrual cycles, which is consistent with the risk analysis 
that shows no elevated risk for younger women using OCs as discussed previously. 
This is also consistent with the pattern analysis by Leinster and Whithouse that shows 
no significant difference in the patterns of pre-menopausal women who have or have 
not used OC. Work by Spicer et al show that when using a reduced estrogen and 
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progesterone OC regimen (reduced to the levels of that found in postmenopausal 
women) there is an associated serial decrease in mammographic density, where the 
authors believe the density mirrors reduced risk. 

[01 29] Other evidence indicates the characteristics of breast tissue are determined early 
in life. Increased placental (the main estrogen producing organ during pregnancy) 
weight may be a perinatal indicator of high risk pattern. This bodes well with work 
that addresses the effects of the proliferation of rudimentary breast epithelium and 
influences due to endogenous hormones and testosterone in fetal development. 

[01 30] The epithelial content decreases with age varying 1 5%-5% by volume from young 
to older women, respectively, but the functional form of this decay is not understood, 
and the relation with what is observed on the radiograph is not clear. Moreover, a 
case-control study of women with a history of benign breast biopsy shows that age, 
menopausal status, and years since last birth influence epithelial percentages, 
whereas oral contraceptive use and body morphology characteristics influence stromal 
proportions; the work also shows decreasing proportion of epithelial content with age. 
This analysis was at the biopsy site and may not apply across the breast but does 
provide circumstantial evidence that the density analysis of mammograms provides a 
crude scenario of the dynamics at play when used in isolation. Likewise, the work 
provided by Russo et al also suggests that densities as measured on mammograms 
are not equal, since the predisposition for BC depends on the lob aggregate, which is 
dependent on the reproductive history. 

[01 31] A critical epithelial-stromal ratio may be important. Likewise, when assessing 

density, there may be an age related correction or conversion function that will allow 
better risk estimates. That is the amount of dense tissue that is at risk may be found 
by applying an age dependent adjustment to the actual image density measurement. 
The evidence indicates the epithelial content decreases with time, but may not be 
resolved in the radiograph since it (the epithelium) is of limited volume to begin with 
and may be masked by stromal content. 

[0132] 

The volume changes during the menstrual cycle noted in the MR work may hold 
some clue. As reported the parameter changes due to MT vary greatly within each 
individual. If in some way this change is more closely related to one of the dense 
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constituents than the other, perhaps a better assessment of the true tissue at risk can 
be made. The research by Bernardes et al shows that tamoxifen has a definite impact 
on reducing cell proliferation for those on tamoxifen therapy as measured over one 
menstrual cycle. The research in MT and mammography indicates the changes in 
mammographic density due to menstrual phase for normal women are observable in 
standard x-ray images for some cases, although quantitative methods were not 
applied for image analysis. These changes are consistent with the MR volumetric work 
discussed above. It is possible that over many menstrual cycles there is a cumulative 
effect due to tamoxifen therapy, or a similar type therapy, which manifests over 
extended periods as a decreased density for some women, in particular for lean 
women. 

[01 33] As pointed out by Hutson et al, mammary carcinoma more often develops in the 
upper outer quadrant of the mammary gland, but the related analysis shows that the 
epithelial content is greater in the lower half. This suggests that there may be a spatial 
aspect of density analysis that has not been considered in the risk analysis as of yet. 
Density may be a risk that has a spatial dependence. Similarly, other work also shows 
a spatial dependency in tumor growth. It may be wise to consider this as a spatial risk 
factor when considering CAD development and explore methods of adjusting 
parameters dependent on spatial location. 

[01 34] Ljtt|e researc h has been aimed at understanding the interplay between diet, and or 
exercise with serial pattern change. Boyd et al showed that women on a reduced fat 
diet showed no significant density or dysplastic change after one year. However, 
subsequent research indicates that the adoption of a reduced fat and increased 
carbohydrate diet influences radiographic appearance. This work, that followed 
women initially with greater than 50% dense breast, shows decreases in breast area, 
primarily accounted for by weight loss, and loss of dense area and percentage area 
attributed to the intervention diet after two years with the effect more pronounced in 
women going through menopause. Subsequent work of similar dietary intervention of 
women going through menopause shows significant reduction in dense area, 
attributable to changes in cholesterol, total fat, and saturated fat intake, and changes 
in percentage density, attributable to changes in cholesterol and saturated fat intake, 
in addition to showing that women with FH of BC show less reduction that those 
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without. Other work shows little or no association with alcohol intake and density 
proportions. This is consistent with alcohol-risk research that shows no increased risk 
for alcohol users less than 50 years of age and increased risk for those over 50. 

[01 3 5] Cross-sectional (case-control) research indicates a positive relation with saturated 
fat (or total fat) intake, and inverse association with fiber and carotenoid intake with 
homogeneous or nodular densities, but no regular association with high risk Wolfe 
pattern. This work also shows no significant relation with pattern and polyunsaturated 
fat or cholesterol intake. Likewise, subsequent work demonstrates a positive 
correlation with fat intake and women with DY patterns. However one report indicates 
increased breast density with decreased intake of saturated fat in pre-menopausal 
women. Other work indicates an association with protein (post menopausal women), 
carbohydrate, and meat (post menopausal women) intake with high risk Wolfe 
patterns and no significant relation with fat or vitamin intake with pattern. In a case- 
control study of ethnicity and mammographic density evidence indicates an inverse 
relation with soy intake and breast density. 

[01 36] Only one case-control study was located that addresses physical activity and 
breast density. When controlling for confounding factors, this work shows a weak 
inverse association with density and moderate physical activity with the strongest 
association found in young women in the highest category of activity. The inference 
could be made that exercise induced morphological changes can influence 
mammographic patterns through alternative avenues such as breast size and body 
style as discussed previously. 

[01 37] In general, Hoffman-Goetz et al provide a good review on the association between 
exercise and BC that includes reproductive, hormonal and endocrine mechanisms. The 
authors indicate that the evidence suggests an inverse relation with exercise and BC. 
This is in agreement with other reviews indicating physical activity reduces the BC risk 
in particular lean, parous and pre-menopausal. 

[01 38] Setting the age-related screening commencement argument aside, the screening 
pool is at the age where significant proportions of the patterns or densities are most 
likely to shift, more often to a lesser degree. For robust automated sequential image 
analysis, these shifts should be understood and incorporated into the analysis 
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scheme. A quantitative measure for "normal" aging associated with density 
measurements or tissue composition analysis obtained from radiographs has not been 
researched to any great extent. 

[01 39] Serial image analysis, in many cases, is synonymous with what is generally termed 
change detection. Change detection has its roots in the (non-medical) image analysis 
of similar scenes acquired at different times. An intuitive method for finding change, 
which is productive in many cases, is sequential image subtraction, which may be 
considered as a time derivative. It follows that large changes result in large signals in 
the time derivative image. Some applications include the study of scenes acquired 
from aerial or satellite surveillance. Techniques such as geometric image warping, 
rigid rotation, and translational alignment (often termed image registration) are 
sometimes employed to correct for atmospheric interference or for differing 
inclination angles. 

[0140] More recently, change detection (serial analysis) applications have found their way 
into the medical imaging analysis arena with applications in DM. it is often inferred 
that geometric warping is an important initial correction that should be applied prior 
to automated analysis methods. This inference is made because it is assumed that 
serial mammographic images will most likely result from differing x-ray projections 
due to plate compression and patient positioning, and is probably a true assumption. 
However, this author does not know of any work that provides a quantitative analysis 
as to the extent of the problem. The work by Kostelec et al also shows that anomalies 
results when applying warping techniques to dissimilar regions assumed the same. 
The constraints or boundary conditions for the warping procedure are that sequences 
of images (a serial pair) are the same; the serial pair only differs due to the 
projection-position distortion. The work presented in this review suggests that this is 
an assumption that is most likely not applicable in many cases. Additional temporal 
effects must also be accounted for in the sequential image analysis that applies to 
pre- and post menopausal women separately. 

[0141] 

For pre-menopausal women, optimal imaging acquisition should include 
menstrual timing, since volumes and density characteristics may change. That is for 
automated analysis, the images should be acquired at the same stage of the 
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menstrual cycle, and if possible during the follicular stage when the breast is less 
dense. Also OC use is a consideration, but the evidence is not clear as to the influence 
on the breast image. This may follow from noting that the circulating hormones due 
to the natural cycles are similar to those during OC use. 

[01 42] For postmenopausal women, the start of HRT therapy should also be considered. 
From a pragmatic stance, the evidence shows that if a change in mammographic 
structure takes place, it does so in the short term. The suggestion here is to update 
the baseline image at about the 24- month mark for women starting this therapy. 
Likewise, the evidence suggests that discontinuing HRT encourages the 
commencement of the normal involution process. 

[01 43] Likewise, there are other factors to consider such as weight change, or change in 
exercise habits, and the general aging in terms of breast density influences. Although 
limited in scope, the research indicates that body morphology is related to breast 
tissue, and that a well designed serial analysis should attempt to track this 
information. 

[0144] There seem to be distinct sources of error that propagate through pattern-risk 
analysis (a) the methodology applied to estimate parenchymal characteristics for 
example Wolfe pattern, BIRADS, or density analysis, (b) the masking effect, (c) possibly 
age related cancer growth rates, (d) the time at which the pattern was assessed 
(present, past or some combination of the past with the present ), (e) the nature of the 
tissue at risk, since dense tissue is a composite of epithelial and stromal content, and 
(f) perhaps there is a spatial dependency for dense tissue and risk. 

[0145] p or tne past decade or so many researchers have developed novel approaches for 
CAD with many methods based on detecting or enhancing the abnormalities. For the 
most part these methods have been applied to digitized film images. The general 
definition of CAD is widened in the present invention to include understanding the 
complicated and entwined factors discussed in this specification. That is, CAD should 
also include developing methods for quantifying the factors discussed here, in 
addition to detection, classification, and teleradiology applications. With the advent of 
real time digital images (the FFDM) system, it may be possible to quantify the 
temporal influences of involution or aging, HRT, environmental, other hormonal 
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interactions, and the pattern-risk assessment and link these factors with the actual 
image measurements. 

[01 46] The FFDM data collection system directly supports the serial analysis of involution 
and time related density change and risk analysis. The automated density analysis 
methods, discussed above, should permit the real time risk assessment, due to dense 
tissue, over extended intervals. This should enable the tracking of large cohorts over 
many years with density measurements made at each screening session for the 
participants. In this fashion, in conjunction with the radiologist's assessment, more 
accurate risk assessments relative to density should result. 

Summary of Invention 

[0147] The present invention is a method of determining breast cancer risk including the 
steps of establishing a risk probably value associated with a patient, the risk 
probability value calculated from an array of risk factors associated with breast cancer, 
applying a computer algorithm adapted to find abnormalities in the patient"s 
mammogram, and increasing the tolerance level for false positive results in the 
computer algorithm responsive to a higher probability value associated with the 
patient and decreasing the tolerance level for false positive results in the computer 
algorithm responsive to a lower probability value associated with the patient. 

[01 48] The risk factors may include relative risk values, odds ratio values or absolute risk 
values. An advantage of the present invention is its utility in automated screening. 
Accordingly, an embodiment of the invention includes the steps of obtaining a 
patient-specific breast tissue density value derived by automated means from the 
patient's mammogram and integrating the breast tissue density value in the risk 
probability value. Mammograms generating a positive result for breast cancer are 
flagged for additional analysis. A recommended course of action may be automatically 
generated wherein more invasive procedures are recommended responsive to the 
higher probability value and less invasive procedures are recommended responsive to 
the lower probability value. 

[0149] 

In another embodiment, a data entry interface adapted to input the array of risk 
factors associated with the patient is provided and the patient's mammogram is 
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digitally acquired. The array of risk factors is stored on an electronic storage medium 
communicatively coupled to the digitally acquired mammogram. The algorithm is 
applied to the mammogram to find abnormalities. Optionally, mammograms 
associated with abnormal risk findings are electronically presented with computer 
aided enhancement. 

[01 50] The array of risk factors includes at least one factor selected from a group of 
factors including age, racial background, geographic background hormonal data, 
breast size, weight and height, pregnancies, breast surgeries, breast water content, 
transverse relaxation time, family medical history, previous biopsies, length of 
reproductive years, menopausal status, parity, age of menarche, age of menopause, 
involution characterization, density time dependency, density dependent texture, 
dietary factors, abnormality spatial location and physical activity. 

Brief Description of Drawings 

[01 51] For a fuller understanding of the nature and objects of the invention, reference 
should be made to the following detailed description, taken in connection with the 
accompanying drawings, in which: 

[0152] FIG. 1 shows Equations 1-6; 

[01 53] FIG. 2 is a diagrammatic view of the method according to the invention; 

[01 54] FIG. 3 is a diagrammatic view of detection criteria damped according to risk, and 

[01 55] FIG. 4 is a diagrammatic view of a plurality of patient history variables that affect 
breast cancer risk as integrated into the present invention. 

Detailed Description 

[0156] 

There are various ways of modeling risk dependent on the experimental design. 
The resulting risk figures of merit are often expressed as an odds ratio (OR), relative 
risk (RR), or absolute risk (AR). In general, RR can be calculated from an incidence 
study (or often termed cohort study), where a group is tracked over time and the 
incidence of the disease is monitored. OR can be calculated from a case-control study, 
which under certain conditions can be used to estimate the RR even though the study 
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design does not strictly support the calculation. The RR estimates from case-control 
studies often agree with the associated RR calculated from cohort studies of the same 
population. The AR can be estimated by understanding the RR associated with a given 
population for a given set of risk factors in addition to understanding the baseline 
incidence for a patient without any of the risk factors considered. 

[01 57] A brief overview of methods used for estimating risk will be presented here. The 
three risk estimates described above can be incorporated into CAD analysis. This is in 
no way a complete discussion of risk analysis, but does provide sufficient background 
to fuse risk analysis with CAD. 

[01 58] Logistic regression techniques can be applied to case-control studies for 

estimating for calculating ORs. This is a useful approach when the observed outcome 
P is dichotomous or restricted to two values representing the occurrence (P = 1) or 
nonoccurrence of some disease (P = 1 ) (P =0) (See Figure 1 , Equation 1 ). The model is 
given by Figure 1 , Equation 2 where the x's are the independent variables and the c's 
are the regression coefficients. Here the independent variables are the risk factors, 
which can be dichotomous or continuous. For example, mammographic density can 
be represented as total density or by proportions or by view; similarly weight and 
height are continuous as well as years of menstruation; family history can be 
represented dichotomously as 1 if present and 0 if not. Standard methods of forward 
or backward selection can be applied to determine the significance of a particular risk 
factor. Or, more generally, conditional logistic regression can be applied to determine 
the model. All standard risk factors are explored, which include age, density, FH, 
reproductive years, parous information, weight, height, and other reproductive 
history, such as age of menarche and age at menopause. It may be the case the two 
models should be developed depending on menopausal status. This work can be 
implemented by standard case-control methods. When sampling a database and 
implementing a case-control study, normally ORs are readily assessable. When certain 
conditions apply, ORs can be used as an approximation for RR. 

[01 59] 

There is a different y for each woman (or patient) that makes up the data set. In 
order to apply this technique, a training data set must be constructed that includes a 
distribution of women resembling that found in the screening environment. The 
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regression coefficients are found by fitting the model with the risk variables measure 
from the training set. The training set includes women with the various risk factors 
with and without the incidence of cancer. Ideally, the cases will be sampled randomly 
from the two pools (those with and without cancer). 

[01 60] One method of estimating the absolute risk of associated with the disease over 

some time interval is to use the Gail model approach. Rather than developing the Gail 
approach specifically, the necessary background from reliability theory is provided, 
from which the many other approaches follow or are derived from. 

[0161] In component reliability analysis, the failure rate X (t) plays an important role. 

There is an analogous expression in survival analysis termed the death rate, which is 
sometimes referred to as the hazard function. Let f(t) represent the probability density 
(pdf) for survival time with corresponding cumulative probability F(t). The death rate 
(or the incidence of some event) is then the rate at which the deaths occur divided by 
the proportion of the surviving population is provided in Figure 1, Equation 3 with the 
obvious substitution. S(t) is referred to as the survivor function (or in reliability theory 
the probability that component will not fail before time t). In the present invention X 
(t) is the rate at which the disease occurs (more appropriately the age) or the age- 
specific incidence rate. The standard solution is provided in Figure 1 , Equation 4. If X 
(t) is under some form of stress (the analogy here is risk), denoted by x, with pdf p(x), 
where x in the most general terms is multivariate, we can then express the observed 
incidence rate in terms of the conditional incidence rate in Figure 1 , Equation 5. The 
proportional hazard model is written in Figure 1 , Equation 6 for a particular set of risk 
factors x, where X ^ (t) is the baseline hazard when all the covariates are zero, and the 
3 's are the regression parameters. Now Equation 6 can be used to derive the 
probability of the incidence of the disease over some time period for a patient with a 
given set of risk parameters. 

[01 62] We have provided for the most part parametric methods for estimating risk. These 
methods are not exclusive. The embodied ideas of this invention are not restricted to 
these methods of calculating risk with a given set of risk factors. The ideas discussed 
here may also be implemented with nonparametric methods. 

[0163] | n Figure 2, a risk probably value associated with a patient is established 10 from 
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the patient"s history 20 and factors derived from the image, such as density, in view 
of weighted risk factors JO , the risk probability value calculated from an array of risk 
factors associated with breast cancer. A computer algorithm 40 adapted to find 
abnormalities in the patient's mammogram 5<?is applied. A standard threshold 60 
represents the tolerance level of the algorithm for identifying abnormalities 
independent of specific patient risk information. The standard tolerance level 60 for 
false positive results is increased in the computer algorithm 40 responsive to a higher 
probability value associated with the patient and decreased for false positive results 
responsive to a lower probability value associated with the patient forming a modified 
threshold 70 for detection. Under the modified threshold 70 , examined tissue is 
identified as normal 80 or abnormal 90 . For results that indicate abnormalities, the 
mammogram 50 is flagged for further examination 100 . 

[01 64] In Figure 3, the modified threshold 70 is variable from zero to one representing 
no risk to high risk respectively. By raising the threshold 70 , false positives are 
decreased. Conversely, by lowering the threshold, false positives are increased. As 
there are finite resources for analyzing mammograms, adjusting the tolerance level 
for false positives optimizes the limited resources available. Patients at high risk 
should have a higher tolerance for false positives while patients at low risk should 
have a relatively lower tolerance for false positives. 

[0165] The risk factors may include relative risk values, odds ratio values or absolute risk 
values depending on the screening situation and patient background. An advantage of 
the present invention is its utility in automated screening. Accordingly, an 
embodiment of the invention includes the steps of obtaining a patient-specific breast 
tissue density value derived by automated means from the patient's mammogram and 
integrating the breast tissue density value in the first in the risk assessment and then 
into the ensuing detection algorithm(s). Mammograms generating a positive result for 
breast cancer are flagged for additional analysis. A recommended course of action 
may be automatically generated wherein more invasive procedures are recommended 
responsive to the higher probability value and less invasive procedures are 
recommended responsive to the lower probability value. 

[0166] 

In another embodiment, a data entry interface adapted to input the array of risk 
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factors associated with the patient is provided and the patient's mammogram is 
digitally acquired. The array of risk factors is stored on an electronic storage medium 
communicatively coupled to the digitally acquired mammogram. The algorithm is 
applied to the mammogram to find abnormalities. Optionally, mammograms 
associated with abnormal risk findings are electronically presented with computer 
aided enhancement. 

[01 67] The array of risk factors shown in Figure 4 includes at least one factor selected 
from a group of factors including age, racial background, geographic background 
hormonal data, breast size, weight and height, pregnancies, breast surgeries, breast 
water content, transverse relaxation time, family medical history, previous biopsies, 
length of reproductive years, menopausal status, parity, age of menarche, age of 
menopause, involution characterization, density time dependency, density dependent 
texture, dietary factors, abnormality spatial location and physical activity. It should be 
noted that the present invention anticipates numerous factors that affect breast cancel 
risk and the enumerated variables should not be considered exhaustive. 

[01 68] The rjsk probability model (either relative or probability of incidence over some 
specified period) can be used to adjust detection thresholds. The cumulative 
probability function or relative risk), or integrated probability distribution function of 
some time interval, has function values (or point value for the latter) ranging from 0 
(no risk) to 1 (high risk). As shown in Figure 3, once acceptable detection sensitivity 
(high true positive detection rates) is achieved without considering risk, the detection 
criteria can be damped according to the risk. Note that for high-risk cases the 
detection sensitivity remains nearly the same (for example, multiply by one), and that 
for low-risk cases the FP rate will be reduced because the sensitivity rate is reduced 
(errors move in tandem up and down). Put simply, for high-risk cases the acceptability 
of FP errors should be increased, and for low-risk cases the FP acceptability should be 
reduced. This should also support normal image detection algorithm performance 
increases. Moreover, since the majority of mammographic screening cases are normal, 
it may be possible to discard a proportion of cases by risk analysis alone. This can be 
implemented in three distinct ways: (1) with OR considerations, (2) through relative 
risk analysis, or (3) through estimating the risk over some specified time interval, 
which leads to an absolute risk figure. All figures of merit have significance, in that 
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the former two may be useful for estimating the relative risk or OR of BC today and 
the latter for estimating the risk in the serial screening case; all figures of merit may 
be used in conjunction with the radiologist's assessment as well as with CAD. 
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